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ABSTRACT 

We present the survey design, data reduction, and spectral fitting pipeline for the VIRUS-P Explo- 
ration of Nearby Galaxies (VENGA). VENGA is an integral field spectroscopic survey, which maps 
the disks of 30 nearby spiral galaxies. Targets span a wide range in Hubble type, star formation activ- 
ity, morphology, and inclination. The VENGA data-cubes have 5.6" FWHM spatial resolution, ~ 5A 
FWHM spectral resolution, sample the 3600A-6800A range, and cover large areas typically sampling 
galaxies out to ~ 0.7i?25- These data-cubes can be used to produce 2D maps of the star formation 
rate, dust extinction, electron density, stellar population parameters, the kinematics and chemical 
abundances of both stars and ionized gas, and other physical quantities derived from the fitting of 
the stellar spectrum and the measurement of nebular emission lines. To exemplify our methods and 
the quality of the data, we present the VENGA data-cube on the face-on Sc galaxy NGC 628 (a.k.a. 
M 74). The VENGA observations of NGC 628 are described, as well as the construction of the data- 
cube, our spectral fitting method, and the fitting of the stellar and ionized gas velocity fields. We also 
propose a new method to measure the inclination of nearly face-on systems based on the matching of 
the stellar and gas rotation curves using asymmetric drift corrections. VENGA will measure relevant 
physical parameters across different environments within these galaxies, allowing a series of studies on 
star formation, structure assembly, stellar populations, chemical evolution, galactic feedback, nuclear 
activity, and the properties of the interstellar medium in massive disk galaxies. 

Subject headings: galaxies: spiral, kinematics and dynamics, techniques: spectroscopic, methods: data 
analysis 



1. INTRODUCTION 

In the context of ACDM cosmology, the forma- 
tion and evolution of galaxies takes place in the 
bottom of the deepest gravitational potential wells 
in the dark matter distribution (dark matter ha- 
los). Gas accretion into these halos and merging 
processes ul timately trigger star fo r mation giving rise 
to galaxies (iBlumenthal et all fl98l iWhite et all 119881 : 
ISomerville fc Primackl Il999f ). Although consensus has 
been reached concerning this picture, the details of 
the baryonic physics behind galaxy formation in the 
centers of dark matter halos are still aggressively de- 
bated. The triggering of star formation an d the vari- 
ables that set the star formation rate (SFR. [Kennicuttl 
19981: iKennicutt fc Evansll2012t iMcKee fc Ostrikerl l2007t 
Lerov et all 120081: iKrumholz et al.ll2009D . the contribu- 
tion from different types of feedb ack processes (AGN, su- 
pernovae, stellar radiation, e.g. iKauffmann et al.lll999t 
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iCroton et all 120061: iThompsonl 12009ft . as well as the im- 
pact of gas a ccretion from the inter-galactic medium 
(IGM, e.g IDekel et all [20091: iFumagalli et all f20TIt 
iDave et ail 120111 : Ivan de Voort et al.ll2011h . at"regulat- 
ing the gaseous budget, structure, chemical composi- 
tion, and kinematics of the ISM, and the role that ma- 
jor and minor mergers as well as secular evolution pro- 
cesses play at shaping galaxie s (|Toomre k, Toomrelll972l ; 
iKormendv k. Kennicuttj[2004T ) . are currently the subjects 
of active research. All these processes play a major role 
in determining how galaxies evolve throughout cosmic 
time, building up their stellar mass and shaping their 
present day structure. 

The detailed manner in which the above physical phe- 
nomena (star formation, gas accretion, feedback, interac- 
tions, and secular evolution) proceed, ultimately deter- 
mines the morphology, kinematics, and chemical struc- 
ture of both the ISM and stellar components of a galaxy. 
We can therefore study these processes by obtaining spa- 
tially resolved measurements of quantities like the SFR, 
the metallicity, the stellar and ionized gas velocity and 
velocity dispersion, the age of the observed stellar pop- 
ulations, the atomic and molecular gas surface density, 
etc., and analyzing how they relate to each other and to 
the global properties of the galaxies under study. In this 
way, we can test current theoretical models describing 
the above phenomena and their impact on galaxy evolu- 
tion. 

Wide field optical integral field unit (IFU) spec- 
troscopy allows the measurement of many of these quan- 
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tities in nearby galaxies. Optical IFU data-cubes, in 
combination with multi- wavelength broad band photom- 
etry, and sub-mm and radio maps of the same galaxies, 
are powerful datasets to study galaxy evolution. 

Integral field spectroscopy of nearby galaxies has been 
somewhat limited in the past, mostly due to the small 
field-of-view of available IFUs. During the last decade, a 
new generation of wide field integral field spectrographs 
like SAURON on t he 4.2m William Herschel Telescope 
(jBacon et al.ll200lD. PPAK o n the 3.5m at Calar Alto 
Observatory (iKelz et alJl2006D. Sparse Pak on the WIYN 
3.5m telescope ifBershadv et"atl 120041 ) . and the Mitchell 
Spectrograph (formerly known and hereafter referred to 
as VIRUS-P) on the 2.7 m Harlan J. Smi th telescope at 
McDonald Observatory (jHill et al.l 120081 ) . have opened 
the path to study nearby systems subtending large an- 
gular diameters on the sky. 

Early surveys of nearby galaxies using IFUs mostly 
focused on studying the kinematics and stellar popula- 
tions of early type systems. These include the S AURON 
Survey (|Bacon et al".ll200lHde Zeeuw et al.ll2002D. and its 
extension, the Atlas3D Survey ( |Cappellari etall l20H 
which by now have mapped hundreds of elliptical and 
lenticular galaxies. Wide field I FU studies o f later type 
disk galaxies include the work of lGanda et al.l (|2006l) who 
used SAURON to observe the central regi ons of 18 nearby 
late- t ype spirals, the Disk Mass Project ([Bershadv et al.l 
12010( 1 which used SparsePak and PPAK to measure Ha 
velocity fields for 146 face-on spirals, and stellar kine- 
matics for a subset of 46 objects, with the aim of con- 
straining the distribution of stellar mass and dark mat- 
ter in disk galax ies, and the PPAK IFS Nea rby Galaxies 
Survey (PINGS. iRosales-Ortega et al.l2010H . which maps 
the disks of 17 nearby disk galaxies. The PPAK IFU is 
currently being used to conduct the Cal ar Alto Legacy 
Integr al Field Area survey (CALIFA, iSanchez et atl 
l2011bl) . a massive project mapping ~ 600 galaxies of 
all Hubble types, selected based on their angular size 
and distance (in order for them to fill the PPAK field- 
of-view). In the future the MANGA SDSS IV projecfl 
will produce optical datacubes for tens of thousands of 
galaxies. A number of IFU studies of galaxies have 
been conducted at high redshift (1 < z < 3), where 
target sizes are well suited to the small fields of view 
of IFUs in 10m class telesco pes (e . g. iGenzel et al.1 
20061: iForster Schreiber et al.l [20061: lLaw et al.l 120071: 
Wright et al.ll2007t iLemoine-Busserolle et al.ll2010D . 

In this work, we present the VIRUS-P Exploration of 
Nearby Galaxies (VENGA), an IFU survey of 30 nearby 
spirals which uses VIRUS-P (currently the IFU with 
the largest field-of-view in the world) to spectroscopi- 
cally map large portions of the disks of these objects. 
The sample spans a wide range in Hubble types, SFRs, 
and morphologies, including galaxies with classical and 
pseudo-bulges as well as barred and unbarred objects. 
Ancillary multi-wavelength data exist for many of the 
targets. This includes HST optical and near-IR imag- 
ing with ACS and NICMOS, Spitzer mid-IR and far- 
IR imaging with IRAC and MIPS, near-UV and far-UV 
imaging from GALEX, and far-IR HERSCHEL data. 
Also CO and HI 21cm maps are available for most of 
the sample. VENGA's potential lies in a combination 

9 http://www.sdss3.org/future/ 



of wide spatial and spectral coverage, good spatial reso- 
lution, and depth. The large 1.7' x 1.7' field-of-view of 
the IFU allows us to typically sample each system out 
to ~ 0.7i?25 by tiling only a few VIRUS-P pointings. 
The size of the VIRUS-P optical fibers (4.24" in diam- 
eter) samples physical scales of ^300 pc at the median 
distance of our targets and makes our observations very 
sensitive to low-surface brightness emission. We obtain 
spectra with a median S/N — 40 per fiber per spectral 
resolution element which allows good measurements of 
stellar absorption and nebular emission line spectral fea- 
tures at the native spatial resolution of the instrument 
over most of the data-cube for every galaxy. 

In comparison to other IFU surveys of spiral galax- 
ies VENGA is factors of a few deeper than surveys like 
PINGS and CALIFA while having a similar wavelength 
coverage and a factor fo two larger spatial resolution. As 
mentioned above a key advantage of the VENGA sample 
is the existance of extensive multi-wavelength ancillary 
data for most of these very nearby galaxies. In particu- 
lar the far-IR, mm, and radio data is essential for study- 
ing the relation between stars and the ISM and the pro- 
cess of star formation. While CALIFA and in the future 
MANGA will produce datacubes for hundreds and tens 
of thousands of spiral galaxies, only few of those targets 
will have such extensive ancillary datasets. In that sense 
the VENGA data can provide a good base to interpret 
results from larger but shallower surveys of more distant 
galaxies. 

In the near future, VENGA will also include a high 
spectral resolution component consisting of IFU mapping 
at R ~ 6000 of the central regions of a su bset of the sam- 
ple u sing the VIRUS-W spectrograph (jFabricius et al.l 
120081 ). The high resolution VIRUS-W observations will 
complement the broad wavelength range VIRUS-P ob- 
servations and will provide a more detailed view of the 
kinematics of gas and stars in the central regions of spiral 
galaxies. This component of the survey will be presented 
in a future publication. In this work we limit ourselves 
to discussing the VIRUS-P data. 

The VENGA data are being used to conduct an exten- 
sive set of studies on star-formation, structure assembly, 
stellar populations, gas and stellar dynamics, chemical 
evolution, ISM structure, and galactic feedback. The 
data will also provide an excellent local universe control 
sample for IFU studies of high-z galaxies. The survey is 
designed with the following science goals: 

• Study the process of star-forma tion on galactic 
scales (jKennicutt &: Evansl 120121 and references 
therein), including the correlations between the 
SFR and the star formation efficiency (SFE) 
with other parameters like gas and stellar surface 
density, metallicity, galaxy dynamics, and stellar 
populations. The ultimate goal is to understand 
what are the relevant parameters setting the SFR 
across different environ ments within galaxies (see 
IBlanc eFld1l2009Ll2012T) . 

• Investigate the assembly of central spheroidal 
stellar components in disk galaxies 
(|Kormendv &: Kennicuttl 12004, and references 
therein). This includes characterizing the dynam- 
ics, stellar populations, and chemical abundances 
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of classical and pseudo-bulges in spiral galaxies, in 
order to constrain their star-formation history and 
understand their origin. The goal is to distinguish 
between different evolutionary paths that might 
give rise to these structures (secular evolution, 
galaxy-galaxy interactions, etc.). 

• Provide detailed observations of radial gas inflow in 
the central parts of disk galaxies, induced by barred 
and oval potentials (|Kormendv fc Kennicuttll2004L 
and references therein). This includes studying the 
velocity field of ionized gas and stars in the regions 
influenced by the presence of bars and ovals, and 
also the effects of induced gas condensation and 
shocks in the local star formation efficiency. 

• Construct two-dimensional maps of the stellar 
and gas phase m etall icity of spiral galaxies (e.g . 
iMoustakas et alJl2010URosales-Ortega et aHl20illf ). 
These maps will allow the study of radial abun- 
dance gradients measured with exquisite detail, the 
dispersion in abundances as a function of galacto- 
centric radius, and deviations from axisymmetry in 
the two dimensional metallicity distribution. Com- 
paring these measurements to chemical evolution 
models will help constrain the chemical enrichment, 
gas accretion, and star-formation history of disk 
galaxies in the local universe. 

• Use nebular emission line diagnostics to study the 
nature of the ionizing sources influ encing differ- 
ent parts of the disks of spirals (e.g . iKewlev et al.l 
12001 lHc3l200fit IHaffner et all [20091 and references 
therein). This includes the study of low luminos- 
ity active galactic nuclei (LLAGN), and their im- 
pact on the physical conditions of gas in the central 
parts of galaxies, and the study of the structure, 
metallicity, and kinematics of the diffuse ionized 
gas (DIG) in the ISM. In particular, for the two 
edge-on systems in the sample, we will also be able 
to constrain the properties of the DIG as a function 
of distance above the mid-plane. All these studies 
will provide insight about different feedback pro- 
cesses at play in star-forming disk galaxies, and the 
transport of UV photons in the multi-phase ISM. 

• Study the distribution of stellar mass and dark 
matter in spiral galaxies, by using a combination 
of the VENGA gas and stellar velocity fields and 
constraints o n the M/L ratio from stellar pop- 
ulatio ns (e.g. iBershadv et al.l l2010t lAdams et al.l 
I2012D . The data should allow us, in principle, to 
set constraints on the shape of the dark matter halo 
density profile on these systems, which we can use 
to test the predictions of ACDM models. Of par- 
ticular importance for these measurements is the 
high spectral resolution VIRUS-W component of 
the survey. 

Throughout this paper we use the VENGA data on 
the face-on SA(s)c galaxy NGC 628 to exemplify our 
methodology. This nearby grand-design spiral galaxy 
has been extensively studied in the past and IFU spec- 
troscopy has been obtained in its central region by 
iGanda eF al. ( 2006) and over a larger area than presented 



here , although to a shallower de p th, by the PINGS sur - 
vey (|Rosales-Ortega et al.l 120111: iSanchez et all 1201 lbj) . 
In this work we limit ourselves to presenting the meth- 
ods used to observe, reduce, and analyze the VENGA 
IFU data, and we postpone the presentation of scien- 
tific results on NGC 628 to separate publications, includ- 
ing an accompanying paper studying the radial profile of 
the CO to H2 conversion factor (Xqo) and its relation 
to the physical co nditions across the disk of the galaxy 
(Bl anc ct al. 2012). We also expect to present and make 
public the full VENGA dataset in a future publication 
once the processing and analysis of the whole survey is 
complete. 

In §2 we describe the survey design, including a de- 
scription of the VENGA sample and the observing strat- 
egy adopted to conduct the survey. We present the 
VIRUS-P observations of NGC 628 in §3. The data pro- 
cessing pipeline and construction of the final VENGA 
data products (i.e. reduced and calibrated spectral data- 
cubes) is presented in §4, followed by a description of the 
techniques used to fit the spectra, measure stellar and gas 
kinematics, and extract emission line fluxes (§5). In §6 
we fit the stellar and ionized gas velocity fields of NGC 
628 using a harmonic decomposition technique and intro- 
duce a new method to measure the inclination of nearly 
face-on systems based on the matching of the stellar and 
gas rotation curves using asymmetric drift corrections. 
Finally, we present our conclusions in §7. Throughout 
the paper we adopt a standard set of ACDM cosmologi- 
cal param eters, H a = 70 km s ~ 1 Mpc~ 1 , f2jvf = 0-3, and 
0\ = 0.7 (punklev et al.ll2009T) . 

2. SURVEY DESIGN 

In this section we present and characterize the VENGA 
sample and we discuss the physical properties of the tar- 
get galaxies. We also describe the observing strategy 
and the instrumental configurations used to execute the 
survey. 

2.1. The VENGA Sample 

Table [1] presents the galaxies observed in VENGA and 
lists their main properties. Targets are chosen to span 
a wide range in Hubble types, from SO to Sd, a wide 
range in inclinations from face-on to edge-on systems, 
and they include both barred and unbarred objects. The 
sample also spans a wide range in right ascension in order 
to allow observations to be carried out throughout the 
whole year, and all objects have declinations 5 > —10° 
to make them accessible from McDonald Observatory. 
Figure Q] presents Digitized Sky Survey (DSST^I) cutouts 
for all the galaxies in VENGA. Overlaid are the 1.7' x 
1.7' VIRUS-P pointings obtained on each galaxy. While 
VENGA is designed to map galaxies out to ~ 0.7i?25, 
for NGC 3198, NGC 4569, NGC 4826, NGC 5055, and 
NGC 7731, only a central pointing was observed due to 
observing time constraints. 

10 The Digitized Sky Surveys were produced at the Space Tele- 
scope Science Institute under U.S. Government grant NAG W- 
2166. The images of these surveys are based on photographic data 
obtained using the Oschin Schmidt Telescope on Palomar Moun- 
tain and the UK Schmidt Telescope. The plates were processed 
into the present compressed digital form with the permission of 
these institutions. 
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Fig, 1. — Digital Sky Surv ey cutouts of the 30 galaxies in the VENGA sample. The targets are ordered by Hubble type (taken from RC3, 
dc Vaucoulcurs ct al. (1991)) from earlier to later. White boxes show the VIRUS-P 1.7' X 1.7' pointings observed on each galaxy. 



Since one of the goals of VENGA is to study the origin 
and properties of stellar spheroids in the inner parts of 
disk galaxies, we included objects showing both classical 
bulges and pseudo-bulges (|Kormendv fc Kennicuttll2004 
and references therein) . To distinguish between these two 
types of stellar structures we adopt a criterion based on 
the Sersic index of t he spheroidal componen t (n^). Fol- 
lowing the results of iFisher fc Drorvl (|2008|) we adopt a 
limit of riB = 2. We consider classical bulges those with 
tib > 2 and pseudo-bulges those with with ub < 2. Ta- 



ble[2]presents the bulge-to-total light fractions (B/T) and 
us values for 18 of the 30 VENGA galaxies taken from 
Dong fc De Robertisl (|2006ft . IFisher fc Drorvl (|2008ft . and 
Weinzirl et al.l ()2009l ) . Of the galaxies for which we found 
bulge Sersic index measurements in the literature, 50% 
are classified as classical bulges and 50% as pseudo-bulges 
using the ub — 2 criterion. 

In order to understand how well the galaxies in our 
sample represent the overall population of galaxies in 
the local universe, we compare their stellar mass (M*) 
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Fig. 2. — Stellar mass versus star formation rate for the VENGA 
galaxies with SFR measurements in Table [3] (red circles) , and star 
forming galaxies in the SDSS MPA/JHU DR7 catalog (black con- 
tours). The red and black histograms show the distributions for 
the VENGA and SDSS galaxies respectively. The stellar mass his- 
togram includes the VENGA targets without SFR measurements. 
NGC 628 is shown as a blue circle. 
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Fig. 3.— Histogram of the logarithm of the VIRUS-P 4.24" fiber 
size in physical units for each galaxy in the VENGA sample (black), 
given the distances adopted in TablelTI Also shown is the size of the 
5.6" FWHM PSF of the final data-cubes (red).The vertical dashed 
lines marks the median native and effective spatial resolutions of 
290 pc and 383 pc. 



and SFR distributions to that of ~ 8 x 10 5 galaxies 
at z < 0.2 from the Sloan Digital Sky Survey (SDSS, 
I York et alJ l 2CjOCh MPA/JHU DR 7FI catalog of star for- 
mation rates|BrinchmanneFan[20M) . For the VENGA 
galaxies, we estimate M* from their i-T-band lu minosity 
by as suming a mass-to-light ratio of Tk = 0.78 ()de Jond 
1996). The if-band luminosities are computed using the 
distances reported in Table [T] and the total if-band ap- 
parent magnitudes from th e 2MASS Large Galaxy Atlas 
(LGA, I.I arret t ct al.1l200l except for NGC 1042, NGC 
3147, NGC 3949, NGC 5981, NGC 7479, and NGC 7331 
which are not included in the LGA so their magnitudes 
are taken from the 2MASS Extended Source Catalog 
(jJarrett et al.l 12000). We do not correct the luminosity 
for dust extinction but we expect this effect to be small 
(~10% in the if-band). 

We were able to find integrated SFR measurements 
in the literature for 24 of the 30 galaxies in th e 
VENGA sample . Thes e are taken from iLee et al.l (|2009D . 
IKennicutt et all (j2003[ ). and lThilker et al.l (|2007D in that 
order of preference when multiple measurements are 
found. The stellar masses and SFRs of the VENGA 
galaxies are reported in Table [3J Figure [5] shows the 
VENGA and SDSS galaxies on the M* versus SFR 
plane. Our sample spans a range in SFR from 0.2 to 39 
Mq yr _1 , and is distributed in this parameter similarly to 
the star forming SDSS galaxies. In term of stellar mass, 
the VENGA galaxies span a range between 8 x 10 8 M Q 
and 3 x 10 n M Q , but with 93% (28/30) of the sample hav- 
ing M* > 10 10 M . By comparing the stellar mass distri- 
bution to that of the SDSS star forming galaxies it is evi- 

11 http:/ /www. mpa-garching.mpg.de/SDSS/ 



dent that our sample is biased towards the high-mass end 
of the local star forming population. Figure [2] shows that 
our targets fall both on and below the M+ — SFR main se- 
quence of star forming galaxies, so starburst systems are 
not present in the sample. The VIRUS-P Investigation 
of the Extreme Environments of Starbursts (VIXENS) 
project is conducting a similar study to VENGA but on 
a sample of 15 interacting/sta rburst galaxies in the lo- 
cal universe (jHeidermanl 120121 ) . VENGA is primarily a 
survey of massive spiral galaxies. Dwarf galaxies, early 
type galaxies, and starburst are not well represented in 
the sample. 

Finally, it is important to characterize the spatial res- 
olution we can achieve with VIRUS-P at the distances 
of our targets. Table [1] reports the physical scale in pc 
corresponding to one arcsecond at the distance of each 
galaxy. The IFU fibers subtend a 4.24" diameter on sky 
which corresponds to a spatial resolution between 80 pc 
for the closest galaxy (NGC 6503) and 1020 pc for the far- 
thest (NGC 5981). As will be described in §4.7 the PSF 
of the final data-cubes is given by the convolution of the 
native instrumental resolution (set by the fiber size) with 
the seeing and a gaussian weighting filter used to com- 
bine different frames at different dither positions. This 
translates into an effective PSF FWHM of 5.6" in the 
final VENGA data-cubes which corresponds to a range 
in spatial resolution of 109 pc - 1.35 kpc, and a median 
resolution of 383 pc. Figure [3] presents a histogram of 
the physical scale in pc of the VIRUS-P native fiber size 
(black) and the final data-cube PSF FWHM (red) for 
our targets. 

2.2. Observing Strategy 
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The VIRUS-P IFU is a square array of 246 optical 
fibers (each 4.24" in diameter) which samples a f .7' x 1.7' 
field-of-view (FOV) with a f/3 filling factor. Three 
dithers provide full coverage of the FOV. The instrument 
images the spectra of the 246 fibers on a 2k x 2k Fairchild 
Instruments charged-coupled device (CCD) with 15 /im 
pixels. The CCD electronics deliver 3.6-4.3 e~ read-noise 
depending on the read-out mode used. Fiber spectra 
have an approximately gaussian spatial profile of ^5 pix- 
els FWHM, and are stacked vertically on the CCD at ap- 
proximately 8 pixels apart from each other, minimizing 
cross-talk between fibers. 

For each VIRUS-P pointing we observe 3 dithers at rel- 
ative positions (Aa, AS) = (0.0", 0.0"), (-3.6", -2.0"), 
and (0.0", —4.0") from the origin to ensure full coverage 
of the FOV. Therefore in each pointing we obtain spectra 
for 738 independent regions. Depending on the angular 
size of the targets we observe up to three pointings on 
each galaxy providing full coverage of the central part of 
the galaxies and a sampling of the disk of most targets 
out to ~ 0.7i?25 (see Figure [T] and last column of Table 
[T|). Overall the VENGA survey consists of 60 individ- 
ual pointings composed of 3 dithers each, amounting for 
spectra of ~ 44, 000 independent regions (typically a few 
100 pc in diameter) across the disks of the 30 galaxies in 
the sample. 

The spectral range of VIRUS-P can be adjusted within 
the 3600-6800 A window. The instrument has a set of 
volume phase holographic gratings which provide differ- 
ent spectral resolutions and wavelength coverages. For 
VENGA we use the lowest resolution grating (VP1) 
which provides a spectral resolution of ~ 5 A FWHM 
and coverage over a ^2200 A wide spectral window . We 
observe each galaxy in a blue setup (3600-5800 A) and 
a red setup (4600-6800 A) , obtaining full spectral cover- 
age in the 3600-6800 A range and doubling the effective 
exposure time in the overlap region (4600-5800 A). All 
the data are taken with lxl binning, which translates 
into a spectral dispersion of ~ 1.11 A pixel -1 , except for 
some early observations of the central pointing of NGC 
5194 whic h were taken using 2x1 binning in the spectral 
direction (|Blanc et al.ll2009f i. 

In terms of depth, the goal of VENGA is to obtain 
spectra that reaches a median S/N ~ 40 in continuum 
per spectral resolution element (FWHM) per fiber across 
each galaxy. For most scientific applications this high 
S/N per fiber requirement prevents us from having to 
bin our data over many spatial resolution elements there- 
fore ensuring good spatial resolution throughout most of 
the maps. Exposure times for different galaxies were 
scaled using their effective .B-band surface brightness 
taken from the RC3 catalog (|de Vaucouleurs et al.lll99~il 
Table [T]) and typically range from 45 min to 3 hr per 
dither. The three dither positions and two wavelength 
setups used imply effective exposure times ranging from 
4.5 hr to 18 hr per pointing. 

3. OBSERVATIONS 

The VENGA survey is still in the phase of data ac- 
quisition. Observations of all targets in the red spectral 
setup started in April 2008 and were completed in July 
2010. Blue setup observations started in September 2010 
and we expect them to be completed during 2013. Table 



0] lists all the observing runs we have conducted as part 
of VENGA, the instrumental setup used, the number of 
nights observed, and the galaxies for which data were 
obtained. 

As mentioned above the typical exposure times for each 
dither range from 45 min to 3 hr, typically divided in 
shorter 15 to 25 min exposures. When conditions are 
not photometric we increase the exposure times to en- 
sure reaching the desired depth. Because of the large 
angular diameter of our targets, during most observa- 
tions the VIRUS-P IFU never samples regions of blank 
sky. Therefore, off source sky exposures are necessary to 
measure and subtract the sky spectrum from the science 
data. We obtain 5 min sky frames bracketing each sci- 
ence exposure. The off-source frames are taken 30' north 
of each galaxy in fields that have been confirmed not to 
contain extended sources. 

Bias frames, arc lamps, and twilight flats are obtained 
at the beginning and end of each night. For the red setup 
we use a combination of Ne+Cd comparison lamps and 
for the blue setup we use Hg+Cd. These combinations 
of lamps provide a good set of strong lines over the full 
spectral range of each setup allowing a good wavelength 
calibration with minimal extrapolation towards the CCD 
edges. 

During most nights we obtain data for one or two 
spectro-photometric standard stars using the six-positio n 
fine dithering pattern presented in iBlanc et al.l (2009) . 
As described in §4.5 and §4.6 standard stars are used to 
perform the relative flux calibration while the absolute 
flux level is calibrated against broad-band images. Dur- 
ing some observing runs the spectra of 1 to 3 out of the 
246 fibers fell off the CCD due to camera and grating 
alignment problems. This translates in a lack of spectra 
for a few fibers in the corner of the field-of-view which 
does not affect the data significantly. 

3.1. NGC 628 Data 

In this paper we present the VENGA data over 3 
VIRUS-P pointings on the face-on Sc galaxy NGC 628. 
A summary of the data are provided in Table Observ- 
ing conditions were variable between different runs and 
within different nights during the same observing run, 
ranging from photometric to partly cloudy conditions 
with average atmospheric transparency down to ^60%. 
The seeing (as measured from the guide star in the co- 
focal guider camera of the instrument) ranged between 
1.3" and 4.0" (FWHM) with a median of 2.0". 

Table presents a summary of the NGC 628 VENGA 
data after rejecting 12 out of 207 frames (6%) which are 
affected by pointing errors or catastrophic sky subtrac- 
tion problems (see §4.4 and §4.6). For each dither in each 
pointing on the galaxy we list the total on-source expo- 
sure time, the number of frames, the average seeing of 
the frames, and the mean atmospheric transparency (as 
measured in §4.6). Overall we obtained 84.2 hours of ex- 
posure on NGC 628 which translates into an extremely 
deep spectral data-cube for this galaxy. The VENGA 
data on NGC 628 is roughly a factor of two deeper than 
for the typical target in VENGA. 

4. DATA REDUCTION AND CALIBRATION 

Data reduction is perfor med using the VACCINE 
pipeline ( Ada ms et al.l 120111 ) in combination with a se- 
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ries of custom built IDL routines. In this section we 
describe the data processing including the fiat-fielding, 
wavelength calibration, and sky subtraction of individ- 
ual science frames, and their flux calibration (both in 
the relative, i.e. across wavelength, and absolute sense). 
We also discuss astrometric corrections applied to the 
data. Both the astrometry and absolute flux calibra- 
tion are based on the comparison of reconstructed broad- 
band images from the VIRUS-P IFU spectra and archival 
broad-band images of the galaxies. Finally we describe 
the method used to combine individual frames taken in 
different wavelength setups and at different positions in 
order to create the final data-cubes and row-stacked- 
spectra (RSS) files for the VENGA galaxies. At each step 
in the data reduction VACCINE and our set of external 
IDL routines propagate a properly calculated error frame 
which is associated with each science frame and used to 
compute the uncertainty in the final spectra. 

Throughout this section we use the NGC 628 VENGA 
data to exemplify the results of different procedures. We 
describe processing steps in the order in which they are 
applied to the data. 

4.1. Basic CCD Processing, Cosmic Ray Rejection, and 
Fiber Tracing 

All individual frames (bias, flats, arcs, sky, and science) 
are overscan subtracted. We combine all the overscan 
subtracted bias files for each observing run (usually ~ 
100) to create an image of any residual bias structure 
which is subtracted from all the flat, arc, sky, and science 
frames. 

We use the LA- Cosmic laplaci a n cos mic ray identifi- 
cation algorithm of Ivan Dokkuml (|200lD to identify and 
mask cosmic rays in the science images. We tune the 
algorithm to be robust enough to identify most cosmic 
rays in the science frames while ensuring that real emis- 
sion lines are not masked. Any residual cosmic rays not 
identified in this pass are removed from the data during 
the combination of multiple frames (§4.6). 

Twilight flats are used to trace the spectrum of each 
fiber on the CCD. VIRUS-P is mounted on a gimbal at- 
tached to the broken Cassegrain focus of the 2.7m tele- 
scope. The gimbal keeps the spectrograph at a con- 
stant gravity vector, making flexure effects on the optical 
path of the instrument negligible. Thanks to this, with 
VIRUS-P there is no need to obtain calibrations at the 
same time and telescope position of the science data. We 
have observed shifts in the positions of fibers on the CCD 
(at the 0.1 pixel level) when large changes in temperature 
occur. VACCINE corrects for these small offsets during 
the flat-fielding stage to properly remove the spatial PSF 
of each fiber from the 2D spectra (§4.3). To minimize the 
magnitude of these temperature induce offsets we use the 
set of twilight-flats which is closest in temperature to the 
average temperature at which the science data are taken 
for both tracing and flat-fielding. Most of the time this 
means the twilight-flats taken at dawn are used. 

We extract the 2D spectrum of each fiber in the sci- 
ence, sky, arc, and flat frames using a 5 pixel aperture 
centered around the pixel containing the centroid of the 
fiber's spatial profile. Using a discrete pixel instead of a 
fractional pixel aperture centered on the trace centroid 
itself avoids having to re-sample the data and conserves 
the noise properties of individual pixels. Using a fixed 



aperture that matches the FWHM of the fiber spatial 
profile implies that we recover 76% of the total flux pass- 
ing through the fiber and makes the effects of neighbor- 
ing fiber cross-talk negligible. We prefer to adopt this 
approach instead of a more sophisticated method to re- 
cover the flux in the wing s of the fiber spatial profile 
(e.g. Gaussian Suppression. lSanchezll2006D as the poten- 
tial ~ 15% improvement in S/N that could be gained in 
the case of perfect deblending and extraction is hampered 
by the systematic flat-fielding uncertainty introduced by 
the correction of temperature induced offsets in the edges 
of the fiber profile. 

At this stage in the reduction VACCINE constructs a 
formal error map which accounts for both read-noise and 
Poisson uncertainty for each pixel in the 2D spectrum of 
each fiber. These maps are properly propagated through- 
out the rest of the reduction (assuming gaussian uncer- 
tainties) and end up providing the weights used when 
combining and collapsing spectra from different frames, 
as well as the flux error in the final data-cube. 

4.2. Wavelength Calibration and Characterization of 
the Instrumental Spectral Resolution 

Arc lamp frames are combined to produce a master arc 
for each night. VACCINE is typically able to match and 
fit ~ 12 emission lines in the blue setup lamps (Hg+Cd) 
and ~ 20 lines in the red setup lamps (Ne+Cd). We fit 
the wavelength solution for each fiber independently us- 
ing a 4th order polynomial. Residuals in the wavelength 
solution show an r.m.s. dispersion of a\ ~ 0.1 A, or a 
tenth of a pixel (~ 6 km s _1 at 5000A). 

We also use the emission lines in the master arc frame 
to characterize the spectral resolution as a function of 
wavelength for each fiber. We measure the FWHM of 
non-blended arc lines by performing single gaussian fits. 
By fitting a second order polynomial to the measured 
FWHM values across the spectral direction for each fiber 
we create a robust map of the instrumental spectral res- 
olution (FWHMj ns ) of each fiber as a function of wave- 
length. Good knowledge of the resolution is essential at 
the time of fitting galaxy spectra with linear combina- 
tions of empirical or synthetic templates which must be 
convolved to the same resolution of the data in order 
to extract meaningful line-of-sight velocity distributions 
(LOSVDs) from the fits. 

The VIRUS-P instrumental resolution is observed to 
change smoothly as a function of position on the detec- 
tor, with values ranging from 4.6A to 6.lA FWHM. No 
variation in the instrumental spectral resolution is ob- 
served between observing runs. In the final data-cube of 
NGC 628 the median instrumental resolution in velocity 
units is o~inst — 123 km s^ 1 . 

4.3. Flat Fielding 

The flat-fielding process in VACCINE is used to divide 
out three different effects from the data: (1) the relative 
fiber-to-fiber throughput, (2) the profile of the fiber PSF 
on the detector across the spatial direction, and (3) the 
CCD pixel-to-pixel variations in quantum efficiency. A 
master twilight flat is created for each night by combin- 
ing the set of flat frames closer in temperature to the 
average throughout the night. First, VACCINE removes 
the signal from the solar spectrum from the master twi- 
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light flat. For each fiber, this is achieved by fitting a 
bspline to the twilight spectra of a set of 60 neighboring 
fibers on the chip (which share a similar spectral reso- 
lution) to create a template solar spectrum (which also 
includes the spectral response of the instrument), and 
then normalizing the spectrum of the fib er of interest 
by th is combined high S/N template (see lAdams et al.l 
1201 If) . Since each fiber provides an independent wave- 
length sampling of the twilight spectrum, combining data 
from a large set of neighboring fibers effectively yields a 
sub-pixel sampled spectrum which after being fit by the 
bspline can be evaluated at the exact wavelength scale of 
the fiber of interest. 

From the resulting normalized flat, VACCINE creates 
a "fiber profile" flat by running a median smoothing ker- 
nel across the spectral direction. This new frame con- 
tains only the relative fiber-to-fiber throughput and the 
fiber PSF spatial profile. Dividing the normalized flat by 
this smoothed version yields a pixel-to-pixel flat that is 
applied to the data. 

As mentioned in §4.1, small (sub-pixel) temperature 
induced offsets in the fiber positions on the CCD are 
sometimes present in the data. This can translate into 
undesired residuals when removing the spatial profile of 
the fibers if the data are divided by a flat-field that is off- 
set in the spatial direction. Residuals can be particularly 
large at the edges of fibers, where the data values are di- 
vided by smaller numbers than at the fiber core. These 
offsets must be corrected for, in order to remove the fiber 
PSF across the spatial direction. To do so, VACCINE 
traces each fiber in the science frame (after running a 30 
pixel boxcar filter across the spectral direction to ensure 
a high S/N measurement of the fiber centroid), and com- 
putes an offset with respect to the same fiber centroid in 
the "fiber profile" flat. These offsets are used to resample 
the smoothed "fiber profile" flat using an optimal sinc- 
interpolation method in order to align it with the data. 
We divide the science and sky frames by this resampled 
smoothed flat frame, therefore removing the fiber spatial 
profile from the data. 

4.4. Sky Subtraction 

Sky spectra are measured by combining information 
from the two off-source sky exposures taken be fore and 
after each science frame. In lBlanc et al.l f|2009T ) we sim- 
ply averaged the before and after off-source frames to 
create the sky frame used for background subtraction. 
While this method worked well on the NGC 5194 data 
presented there, those observations were taken far from 
twilight and under very stable and dark conditions. We 
find that when conditions are not optimal (e.g. close 
to twilight, when clouds are present, or near moon-rise 
or moon-set) the sky brightness can change non-linearly 
with time, making the simple averaging of bracketing 
sky frames insufficient to produce an adequate sky sub- 
traction. We adopt a more sophisticated method to es- 
timate and subtract the sky spectrum from our data, 
which makes use of the temporal information we can ex- 
tract regarding the variability of the sky brightness as a 
function of wavelength, from all the sky frames obtained 
throughout a single night. 

As an example, Figure 0] shows the raw (i.e. not flux 
calibrated) sky spectrum measured from 13 off source 
sky frames taken during the night of November 7th 2008 
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Fig. 4. — Sky spectrum in raw units (before flux calibration) at 
different UT times (color coded) during the night of November 7th 
2008 for the VIRUS-P red setup. 



using the red VIRUS-P setup. The spectra are color 
coded by UT time, with purple at the beginning of the 
night and red at the end of the night. We adopt the 
UT time corresponding to the middle of each exposure. 
The elevated brightness and blue color of the spectrum 
at the beginning of the night is due to the first quar- 
ter moon. Moonset at McDonald Observatory on that 
date occurred at 8.47 hr UT and can be clearly seen as 
a sharp drop in the sky brightness, particularly at blue 
wavelengths. For this night, the darkest skies occurred 
between 9 hr and 10 hr UT. These arc followed by a 
monotonous increase in sky brightness which is steeper 
at redder wavelengths. This brightening and reddening 
is due to both the approachment of twilight and the fact 
that the observations were being done at increasing air- 
mass. Our sky subtraction method uses all this informa- 
tion regarding the wavelength and time dependance of 
the sky brightness throughout each night to correct the 
two bracketing off-source sky exposures so they match 
the sky spectrum at the time the science frame is taken. 

To trace these changes we average the spectra of all 
fibers in each individual sky frame and then divide this 
average sky spectrum into 500 wavelength bins (each cor- 
responding roughly to a spectral resolution element). On 
a night-by-night basis, we use the time stamps of all sky 
frames to construct a nightly "light curve" of the night 
sky surface brightness at the central wavelength of each 
of these 500 bins. Figure [5] shows the relative change in 
sky brightness across the same night presented in Figure 
4, for three different wavelength bins near the blue end, 
middle, and red end of the spectrum (blue, green, and red 
respectively). Measurements for individual sky frames 
are shown as filled circles. For reference, the black open 
diamonds show the change in bolometric brightness (i.e. 
in the total flux integrated across the whole spectrum). 
The trends described in the last paragraph are clearly 
seen. At the beginning of the night the sky brightness 
falls more steeply in the blue than in the red as the moon 
sets, and rises faster in the red than in the blue as we ap- 
proach twilight at the end of the night. Vertical dashed 
lines mark the beginning and end of observations of the 
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Fig. 5. — Relative sky brightness as a function of UT time for 
the same night shown in Figure [4] at three different wavelengths 
(blue, green, and red). Filled circles correspond to measurements 
of the sky brightness from the off-source background frames. The 
beginning and end of observations of the same target are shown as 
vertical clashed lines. Solid color curves show cubic spline fits to 
the sky brightness. The black open diamonds and black solid curve 
show the relative sky brightness averaged over the full spectrum. 

same galaxy. The discontinuity in the sky brightness at 
these times is expected since the telescope is pointed in 
a different direction. We refer to each of these sections 
of the night as "observing blocks" . 

Within each individual observing block, we fit the tem- 
poral evolution of the sky brightness in each wavelength 
bin using a cubic spline (color solid curves in Figure [5]) . 
The best-fit splines allow us to evaluate the overall spec- 
tral shape of the sky emission at any arbitrary UT time 
within an observing block. This allows us to calculate a 
wavelength dependent normalization curve by which the 
two bracketing sky frames must be multiplied in order to 
reconstruct the sky spectrum at the exact UT time of the 
science exposure. Since the normalization factors com- 
puted in the 500 wavelength bins are inherently noisy, we 
fit them as a function of wavelength using a fifth order 
polynomial. This allows us to multiply the sky frames by 
a smooth function of wavelength containing the correc- 
tion without introducing further noise into the observed 
sky spectrum. We apply this procedure to correct the 
before and after sky frames to the UT time of each sci- 
ence frame, and then average the two corrected frames 
to create a single sky background file associated with 
each science exposure. Sky subtraction is performed by 
VACCINE using this composite corrected sky frame. 

Once this final sky frame is available, the method used 
by VACCINE to subtract the sky spectrum is analog to 
the bspline algorithm used to remove the solar spectrum 
from the twilight flats (see §4.3). Briefly, for each fiber in 
the science exposure the sky spectra of a set of neighbor- 
ing fibers in the background frame is simultaneously fit 
using a bspline and is then subtracted from the science 
fiber spectrum after scaling for the difference in expo- 
sure time between science and sky frames. As mentioned 
above, this procedure greatly benefits from the sub-pixel 
sampling of the sky obtained by combining information 
from different fibers having slightly different wavelength 
solutions. 

The quality of the sky subtraction in VENGA is ex- 
cellent. We typically see sky subtraction residuals that 
are fully consistent with Poisson plus read-noise uncer- 



tainties. Larger systematic residuals usually appear at 
the wavelengths of the 4 brightest sky emission lines in 
our wavelength range. These resiudals are caused by the 
fast time variability of these spectral features, which is 
independent of the variability of the sky continuum and 
broad emission bands which dominate our corrections. 
These bright sky lines are masked during the analysis 
of the science data. Rarely (less than 4% of the frames 
in the case of NGC 628), when observations are taken 
under extremely bad observing conditions (usually com- 
binations of clouds and moon, or clouds and proximity 
to twilight), obvious residuals in the sky subtraction can 
be observed in the sky-subtracted science exposures. We 
reject these frames from the dataset. 

4.5. Spectrophotometry Flux Calibration 

We use observations of spec troph ptom e tric s tandard 
stars from iMassev et all (|1988[ ) and [00 (|1990j ) to cal- 
ibrate the VENGA spectra. The method used to con- 
struct sensitivity curves from the IFU ob s ervat ions of 
standard stars is described in I Blanc et al.l ((2009) . The 
only difference with the approach taken in that previous 
work is that here standard stars are only used to perform 
the relative flux calibration as a function of wavelength. 
The absolute scale of the calibration comes from a com- 
parison to broad-band optical images of the VENGA 
galaxies. The absolute calibration is described in the 
next section. 

We calculate sensitivity curves for all standard stars 
taken during each observing run. These curves are nor- 
malized to a common scale and averaged to create a mas- 
ter sensitivity curve for each run. All science frames ob- 
tained during each observing run are multiplied by these 
curves to convert them to units of erg s -1 cm -2 A^ 1 . 
Since the atmospheric transparency at the time of the sci- 
ence and standard star observations is unknown we need 
the absolute scaling described in the following section. 
The error maps for each frame are scaled accordingly to 
the data. By looking at the dispersion between different 
sensitivity curves within each month, we estimate a typ- 
ical uncertainty in the relative flux calibration of ~ 8%. 

4.6. Astrometry and Absolute Flux Calibration 

The pointing of VIRUS-P is done using an offset guider 
camera which images a 4.5' x 4.5' field ~ 9' north of the 
IFU science field. An astrometric calibration between the 
two fields allows the observer to point the IFU by putting 
a guide st ar at specific pixel coordinates on the guider 
detector. lAdams et al.l (|2011l ) found systematic offsets 
of the order of 1" between the same pointings made over 
different observing runs. These offsets are also present 
in our data and must be taken into account before com- 
bining multiple frames. In order to accurately recover 
the astrometry of the VENGA science observations we 
register reconstructed broad-band images of the galaxies 
made from the IFU data to archival broad-band images. 
These comparisons are also used to calibrate the spectra 
in terms of absolute flux. For the NGC 628 blue and red 
IFU data we use the SDSS-II0 DR8 .g-band and r-band 
mosaic^. 

12 SDSS-III is managed by the Astrophysical Research Consor- 
tium for the Participating Institutions of the SDSS-III Collab- 
oration including the University of Arizona, the Brazilian Par- 
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Fig. 6. — Attempted and actual relative positions of the three sets of dithered exposures for the three pointings obtained on NGC 628. 
Stars and solid circles mark the attempted fiducial positions for dithers 1, 2, and 3 (red, green, and blue respectively). Crosses mark the 
actual position at which each exposure was obtained. Open squares and dashed color circles show the average fiber position of the actual 
observations. 



For every science exposure we integrate the spectrum 
of each fiber over the corresponding SDSS transmission 
curve (g or r for blue and red setup data respectively) in 
order to measure the monochromatic flux at the effective 
wavelength of the broad-band filter. Simultaneously we 
convolve the SDSS image of the galaxy with a gaussian 
kernel to match the PSF to the seeing under which the 
VIRUS-P data was taken and we perform aperture pho- 
tometry at the fiducial position of each fiber. We use 
circular apertures that match the fiber size. After cor- 
recting the VIRUS-P fluxes for atmospheric extinction at 
the airmass of the observations we compare them to the 
SDSS fluxes and fit them using the following expression: 

fsoss = Axf VP + B (1) 

where A is a normalization factor recovering the absolute 
flux scale and B recovers any residual background left 
from the sky subtraction process. A perfect background 
subtraction in both the SDSS image and the VIRUS-P 
spectra should translate in B = 0. 

We perturb the fiducial astrometry of the VIRUS- 
P pointing by applying offsets in both right ascension 
and declination in order to minimize the \ 2 °f the fit. 
This registering process provides corrected astrometry 
for each science frame. Given the wealth of spatial in- 
formation encoded in the relative brightness of hundreds 
of VIRUS-P fibers, the registering is very accurate and 
has a typical uncertainty of ~ 0.1". The value of A 
at the registered position provides the absolute scale for 
the flux calibration and the science frame is multiplied 
by this value. This ties our flux calibration to that of 
SDSS which has a zero-point uncertainty of 2% in the 
bands used here. For the 195 individual frames used to 
construct the NGC 628 data-cube, we measure a mean 

ticipation Group, Brookhaven National Laboratory, University of 
Cambridge, Carnegie Mellon University, University of Florida, the 
French Participation Group, the German Participation Group, 
Harvard University, the Instituto de Astrofisica de Canarias, the 
Michigan State/Notre Dame/JINA Participation Group, Johns 
Hopkins University, Lawrence Berkeley National Laboratory, Max 
Planck Institute for Astrophysics, New Mexico State University, 
New York University, Ohio State University, Pennsylvania State 
University, University of Portsmouth, Princeton University, the 
Spanish Participation Group, University of Tokyo, University of 
Utah, Vanderbilt University, University of Virginia, University of 
Washington, and Yale University. 
13 http://data.sdss3.org/mosaics 



(B) = — 2 x 10 18 erg s 1 cm 2 A 1 with a standard 
deviation of ct_b = 4 x 10 -18 erg s _1 cm ~ 2 A" 1 . This 
level of sky subtraction residuals correspond to less than 
2% of the median continuum level in the data and this 
is only an upper limit since the residuals have some con- 
tribution from the error in the SDSS sky subtraction. 

Figure [6] shows the astrometric offsets we measure for 
the three pointings on NGC 628 with respect to the fidu- 
cial dithering pattern. For each pointing the fiducial po- 
sitions of a fiber on dithers 1, 2, and 3 are marked by 
black stars and solid red, green, and blue circles respec- 
tively. Color crosses mark the actual positions at which 
independent science exposures were obtained. These po- 
sitions are measured using the registration method de- 
scribed above. The black squares and dashed color cir- 
cles show the average position for all exposures in each 
dither. 

Overall the VIRUS-P acquisition accuracy is good 
compared to the fiber size. We observe systematic offsets 
from the fiducial dithering pattern in the range 0.1"-4.0", 
with a mean of 1.7". Th is is in good a g reeme nt with the 
1.8" accuracy found by lAdams et al.l (|201 ID ■ As men- 
tioned above, once the VENGA frames are registered to 
the broad-band images and these offsets are corrected we 
estimate an astrometric uncertainty of 0.1". 

When registering the data and building the final 
VENGA data-cubes (next section) we ignore the effects 
of atmospheric differential refraction (ADR). At a typi- 
cal observation airmass of x = 1.2 we expect less than 1" 
of ADR from the blue to the red end of our data. Given 
the 4.2" fiber size of VIRUS-P and the 5.6" FWHM PSF 
of the final data-cubes (§4.7), we decide not to apply an 
ADR correction. 

To evaluate the quality of our spectrophotometric cal- 
ibration we compare the VENGA spectrum of NGC 628 
to the 20" x 20" long-slit d riftsca n spectrum obtained 
by iMoustakas fc Kennicuttl 1)20061 ) as part of the ancil- 
lary data for the Spitzer Infrared N earby Galaxies Sur- 
vey (SINGS, IKennicutt et al.l 120031) . To do so, we use 
the final combined VENGA data-cube described in the 
following section. We integrate the IFU spectra over the 
same 20" x 20" region sampled by the SINGS long-slit 
driftscan. 

The comparison is shown in top panel Figure [7] We ob- 
serve an 18% offset between the two spectra. This offset 
is produced by contamination from the wings of a bright 
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Fig. 7. — Top panel: Comparison between the SINGS 
driftscan spectrum of the cen tral 20" x 20" of NGC 628 (red, 
IMoustakas fc Kennicuttl l20Qfj) and the VENGA IFU spectrum 
(dotted blue) integrated over the same region. The solid blue 
curve shows the VENGA spectrum after masking a foreground 
star which contaminates the measurement. Also shown are the 
residuals between the SINGS and VENGA spectra. Horizontal 
dashed lines mark a ±10% deviation. Bottom panel: Comparison 
between the monochromatic broad-band fluxes measured from the 
integrated VENGA spectrum over the whole datacube (squares) 
and the SDSS g and r and SINGS B and V band images (green, 
orange, light blue and yellow circles respectively). Error bars indi- 
cate the 8% uncertainty in the relative flux calibration. Also shown 
for reference are the integrated VENGA spectrum (red) and the 
broad- band transmission curves (dashed lines). 



foreground star which falls slightly outside the SINGS 
driftscan box, but which significantly contributes to the 
flux inside the box at the 5.6" resolution of the VENGA 
data-cube. Masking the contaminated spaxels (see §6.1) 
brings the VENGA and SINGS spectra to agreement. 
Except for regions on top of strong emission and ab- 



sorption lines, where residuals arising from differences in 
spectral resolution dominate, both spectra show less than 
5% deviations from each other over 80% of the overlap 
spectral range, rising towards the 10% level at the blue 
end of the spectrum. This is well within the 10%-30% 
uncertainty in the absolute calibration and the 3%-4% 
uncertainty in relative fl ux calibration of the SINGS spec - 
trum as estimated in Mousta kas fc Kennicuttl ([2006) , 
and the 8% relative flux error in the VENGA data (de- 
rived from the dispersion between different sensitivity 
curves in §4.5). 

To further test the quality of the VENGA flux calibra- 
tion we compute broad-band monochromatic fluxes from 
the integrated spectrum over the full VENGA datacube 
and compare them to the flux measured over the same 
area in both the SDSS g and r-band images, and the 
SINGS B and F-band optical ancillary image S The 
SINGS images are independent from our flux calibration. 
As shown in the bottom panel of Figure [7] the VENGA 
data agrees with the flux measured in all four images 
within the quoted uncertainties. From these tests we 
conclude that the accuracy of the VENGA absolute and 
relative flux calibrations is good. 

4.7. Final Combined Data-cubes and RSS Files 

At this stage in the data reduction process we have a 
reduced, sky subtracted, wavelength and flux calibrated 
2D spectrum of each fiber (and its associated 2D error 
spectrum) in every individual science exposure. By 2D 
spectrum we mean a non-collapsed spectrum of the fiber 
PSF on the detector, which contains no spatial informa- 
tion about the source, given the fact that optical fibers 
scramble the input light. Also, the position of each fiber 
on the sky is known to 0.1" accuracy thanks to the reg- 
istering process described in the last section. 

Here we describe the methods used to extract and com- 
bine data from different exposures into a final data-cube 
containing a ID spectrum at each position on the sky. 
We also describe the construction of our final row-stacked 
spectra (RSS3) files. These multi-extension FITS3 files 
are described below and contain the spectrum and error 
spectrum, as well as the J2000 equatorial coordinates for 
each spaxel in the data-cube. Also contained in the RSS 
file is information regarding the central wavelength and 
spectral instrumental resolution for every pixel in the 
spectra. 

Data from different nights and observing runs have in- 
dependent wavelength solutions. Therefore, before com- 
bining we need to re-sample all the spectra to a common 
wavelength grid. In VENGA we produce two versions of 
the same data-cube, one with a regularly spaced linear 
sampling of wavelength and another with a regular loga- 
rithmic sampling (i.e. spaced regularly in velocity space). 
The linear version has pixels spaced by 1.1 A (similar to 
the average dispersion in the original data), while the 
logarithmic data-cubes have pixels that are spaced by 
(AA/A)c = 60 km s . We work with spectra in flux 

14 http://irsa.ipac.caltech.edu/data/SPITZER/SINGS/ 

15 RSS files, in which every row contains the spectrum of a 
spaxel, should be distinguished from data-cubes, as they are two 
dimensional instead of three dimensional representations of IFU 
data. They must be accompanied by a position table for all spax- 
els in order to allow for a geometrical reconstruction of the data. 

1B Flexible Image Transport System 
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density units so total flux is properly conserved when 
re-sampling. 

The reason behind producing two versions of the data 
is that while most users will be interested in using the 
linear version for many applications, the spectral fitting 
software used in the following section to extract stellar 
and gas kinematics, as well as emission line fluxes, re- 
quires input spectra that is regularly sampled in velocity 
space. Instead of interpolating the spectra to a linear grid 
for the effects of combining, and later re- interpolating the 
combined spectra to a logarithmic scale, we do both re- 
samplings directly from the original data. In this way 
we avoid the effects of S/N degradation associated with 
extra interpolations. 

For each fiber in each science frame, after re-sampling 
the 2D spectra we have 5 pixels at any given wavelength 
which provide a measurement (with its associated error) 
of the flux density at a certain position on the sky. We 
then define a regularly spaced spatial grid of spaxels. The 
grid has a spaxel scale of 2" chosen to roughly Nyquist 
sample the final data-cube PSF. For every spaxel in the 
grid we compute the final spectrum by combining the 
spectra of surrounding fibers using a Gaussian spatial 
filter and adopting an inverse variance weighting scheme. 
Therefore, the spectrum of each spaxel in the data-cube 
is given by 



/a(A) 



E ^,i,*(A)/r' fc (A) 



(2) 



where /M' (A) is the flux density of pixel k (k = 1 — 5) 
at wavelength A of fiber j in frame i, and the weights 
Wi,j,k are given by 



' J i,j,k 



(A) 



1 



(o-ij.feCA)) 2 



x e 



(3) 



where 0i,j,fc(A) is the error in the flux, dij is the dis- 
tance in arcseconds between the fiber center and the 
spaxel in question, and cr g is the width of the spatial 
Gaussian filter used to combine the data. We use a value 
of <r g = 4.24/2.355 which matches the fiber size to the 
FWHM of the Gaussian filter. The sum in Equation 2 
is performed over all fibers within a maximum distance 
of the spaxel of interest. This distance cut is given by a 
99% drop in the Gaussian filter weight. During the com- 
bination process we apply 3<r clipping rejection to remove 
any cosmic rays that were not masked by LA-Cosmic in 
§4.1. The final spatial PSF of the data-cube is given by 
the convolution of the top-hat 4.24" diameter fiber profile 
with the seeing (2.0" FWHM for the NGC 628 data) and 
the 4.24" FWHM Gaussian spatial filter used to combine 
the data. Figure [5] presents these three components and 
the product of their convolution (red solid line) . The fi- 
nal PSF is well described by a 5.6" FWHM Gaussian and 
is largely independent of the seeing. Direct measurement 
of the FWHM of bright foreground stars in the NGC 628 
data-cube empirically confirms the value of 5.6". 

As mentioned in §4.2, the instrumental spectral res- 
olution at different wavelengths, for each fiber, is ex- 
tracted from the arc lamps. To create a spectral res- 
olution map for the data-cube, we combine the master 
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Fig. 8.— The final PSF of the VENGA data-cubes given by the 
convolution of the fiber profile (solid black line) with the seeing 
(dotted black line) and the Gaussian spatial filter used to combine 
the individual exposures (dashed black line) is shown in red. The 
PSF is well described by a Gaussian PSF with 5.6" FWHM (dashed 
blue line) , and is largely independent of the seeing. 



arc lamp frames associated with each individual science 
frame in the same way as the science data (i.e. using the 
same weights). We then use the method described in §4.2 
to create a map of the instrumental spectral resolution 
from this combined arc. 

The final VENGA data products for each galaxy are a 
data-cube (used mostly for visualization purposes) and 
a RSS file (used for spectral fitting and analysis). The 
data-cube consists of a multi-extension fits file, in which 
each extension corresponds to an image of the galaxy at a 
given wavelength (sampled at either linear or logarithmic 
steps as described above). The data-cubes are in units 
of flux density (erg s _1 cm~ 2 A -1 .). 

The corresponding linear and logarithmically sampled 
VENGA RSS files are stored as multi-extension FITS 
files, and contain the following information in their dif- 
ferent extensions: 

1. Flux density spectrum for each spaxel in units of 
erg s _1 cm~ 2 A -1 . 

2. Error spectrum for each spaxel in units of erg s" 1 
cur 2 A -1 . 

3. Central wavelength of all pixels in the spectrum of 
each spaxel in units of A. 

4. Right ascension and declination of each spaxel in 
units of decimal degrees. 

5. Instrumental spectral resolution (FWHM) as a 
function of wavelength for each spaxel in units of 

A. 

The final VENGA data-cube of NGC 628 covers an 
area of 5.2' x 1.7' and samples the disk of this galaxy 
out to 8 kpc (0.6i?25)- Figure [9] shows an r-band image 



The VIRUS-P Exploration of Nearby Galaxies (VENGA) 13 



o 




arcsec 



Fig. 9. — Map of the r-band flux reconstructed from the VENGA spectral data-cube of NGC 628. Black contours show isophotes in the 
SDSS r-band image of the galaxy (PSF matched to the 5.6" FWHM PSF of the VENGA data which is shown by the white circle in the 
top right corner). Black dots mark the position of each spaxel in the data-cube. This and all maps presented in this work were constructed 
using the PLOT_VELFIELD IDL routine written by Michele Cappellari ( jhttp : // www-astro.physics . ox. ac . uk/ ~mxc /idl/ \ . 



of NGC 628 reconstructed from the final VENGA data- 
cube. All maps presented in this paper have north point- 
ing up and east pointing left. Black dots mark the posi- 
tion of each spaxel in the data-cube and the white circle 
on the top right corner shows the PSF FWHM. The final 
combined spectra was integrated over the SDSS r-band 
transmission curve to create this map. For comparison, 
isophotal contours from the SDSS r-band image (with 
the PSF degraded to 5.6") are overlaid in black. The 
similarity between the two surface brightness distribu- 
tions confirms that the flux calibration and astrometric 
correction of the individual science frames has been done 
properly. The spectra in the final data-cube of NGC 628 
has a median S/N = 117 in the continuum per spectral 
resolution element (FWHM). In the central parts of the 
galaxy we typically have S/N > 300, while the spectra 
in the faintest regions has S/N ~ 12. In terms of surface 
brightness depth, the NGC 628 data reaches a 5a detec- 
tion limit for continuum at a V^-band surface brightness 
of 25.0 mag arcsec -2 (AB). We remind the reader that 
the data on NGC 628 is about a factor of two deeper 
than for the average VENGA target. 

5. SPECTRAL ANALYSIS PIPELINE 

In order to extract emission line fluxes, gas and stel- 
lar kinematics, and information about the stellar pop- 
ulations present in different parts of the galaxies, we 
fit the VENGA spectra using a linear combination of 
templates convolved with a LOSVD plus a set of Gaus- 
sian em ission line profiles. To do the fitting we use the 
pPXF (ICappellari fe Emselleml [200l and GANDALF 
(jSarzi et al.l 120061) IDL routines developed for this pur- 
pose by the SAURON team. In this section we describe 
the fitting process and present an example fit in the NGC 
628 data-cube. We also discuss how the uncertainty in 
the fitted parameters is calculated. In this paper we fit 
the stellar continuum using empirical stellar templates 
(§5.1). The same type of fitting can be performed us- 
ing stellar population synthesis (SPS) templates to ex- 
tract information about the stellar populations present 
in different regions within the galaxies. The SPS fit- 
ting of the VENGA galaxies will be presented in a fu- 



ture publication. Before fitting the spectra, the data is 
corrected for Galactic extinction adopti ng a Milky Way 
(MW) extinction law as parametriz ed bv | Peil (11992t). and 
extinction values from the maps of iSchlegel et al.l |l998[ ) 
(E(B — V)mw — 0.07 mag or Ay mw = 0.23 mag for 
NGC 628). 

5.1. Stellar Kinematics 

We mask the spectrum of each spaxel around regions 
affected by sky subtraction residuals due to bright sky 
lines, regions potentially affected by the nebular emission 
lines listed in Table HI and the edges of each of the indi- 
vidual instrumental setups (red and blue). We then fit 
for the stellar line-of-sight velocity (v*) and velocity dis- 
persion (a* ) with the pPXF soft ware, which uses the "pe- 
naliz ed pixel" fitting technique ( Cappell ari fe Emselleml 
l2004f ) to fit the spectrum as a linear combination of tem- 
plates convolved with an LOSVD. Although the software 
uses a Gauss-Hermite polynomial LOSVD, and allows for 
the fitting of high order terms (h^jhi), for the data pre- 
sented here we only fit for the first two moments v* and 

The logarithmically sampled RSS files are used as in- 
put for pPXF, which requires the spectra to be reg- 
ularly sampled in velocit y space. We use the MILES 
stellar library version 9.1 (ISanchez-Blazquez et al.ll2006t 
iFalcon-Barroso et al.ll2011[) as a source of empirical tem- 
plates. A subset of 48 stars spanning a wide range 
in spectral types (O through M), luminosity classes (I 
through V), and metallicities (—2 <[Fc/H]< 1.5) is used. 
Also included in the subset are horizontal branch and 
asymptotic giant branch (AGB) stars. 

Before fitting, the templates are re-sampled to the 
wavelength scale of the data. We match the spectral 
resolution of all spaxels in the VENGA data and the tem- 
plates to the worst instrumental resolution at any given 
wavelength in the data-cube by convolving the spectra 
with a running Gaussian kernel with a wavelength de- 
pendent width. This translates in a spectral resolution 
of ai ns = 110 km s -1 in the red edge and ai ns — 200 
km s -1 in the blue end of the spectrum. The median 
instrumental resolution across the convolved data-cube 
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is a i ns = 150 km s . To convolve the templates we as- 
sume the corrected MILES library intrinsic resolution of 
2.54A (|Beifiori et al.ll20lTI ). 

Close inspection of fits to very high signal-to-noise 
spaxels in the data-cube (S/N > 100) shows systematic 
residuals caused by template mismatch at the 1% level 
m the 5000A-6800A region, and rising up to 2% towards 
the blue end of the wavelength range. It is important to 
include this systematic uncertainty floor during the fit- 
ting process. Particularly for high S/N spectra where it 
becomes the dominant source of discrepancy between the 
best fit models and the data. We assume a systematic 
uncertainty of 2% at 3600 A, linearly decreasing towards 
1% at 5000 Aand staying flat at this level redwards. We 
add this uncertainty in quadrature to the photometric 
errors before fitting the spectra. 

In order to account for the effect of dust extinction on 
the shape of the continuum as well as systematic differ- 
ences in the flux calibration of the data and the tem- 
plates, during the minimization we fit both additive and 
multiplicative low order Legendre polynomials. The low 
order of the polynomials prevents them from introducing 
features on small scales of the order of the instrumental 
resolution. The polynomials only match the large-scale 
(~ 100A) shape of the continuum in the linear combina- 
tion of stellar templates and the data and they do not 
affect the fitting of individual spectral features. 

We fit the spectrum of each spaxel individually and 
store the kinematic parameters (u*, cr*). We limit the 
wavelength range to the Mg I b absorption feature (5100- 
5300 A) for the purpose of extracting the stellar velocity 
dispersion. We do this to avoid a systematic overesti- 
mation of cr* caused by the broad wings of Balmer ab- 
sorption lines which are subject to rotation and pressure 
broadening. As discussed below, this effect is partic- 
ularly significant in the spectra of star-forming regions 
which have a large contribution from A type stars. We 
also experimented with using other strong absorption 
features like Ca H+K and the G-band, finding that Mg I 
b provides a better agreement with the published high 
precision measureme nt of the velocity disp ersion pro- 
file of NGC 628 by IHerrmann etHI (f2008h (see §6.4) 
while the other features always produced systematically 
higher values of cr*. Issues regarding the overestima- 
tion of cr* when using Balmer lines, Ca H+K and the 
G-b and have been extensively discussed in the literature 
Ce.g.lKormendv fc Illingworthll982l;lBernardi et al.l2003t 
iGreene fc Holl2006f) . The instrumental resolution in the 
convolved spectra at the Mg I b feature is <Ti ns — 145 km 
s- 1 . 

We also perform a second fit using the full wavelength 
range to estimate u*. Our simulations show that while 
the best fit cr* is biased when using the full spectrum, 
the best fit i>* is not and the uncertainty in its value is 
reduced significantly by using the full wavelength range. 
After fitting, we keep the LOSVD parameters fixed over 
the next fitting iteration (§5.3), in which we also fit for 
the emission lines in the spectrum. Kinematic maps of 
NGC 628 are presented in §6. 

5.2. Error Estimation for LOSVD Parameters 

To estimate the uncertainty in the best-fit LOSVD pa- 
rameters (v* and cr*) we fit a series of Monte Carlo re- 



alizations of simulated spectra with different stellar ve- 
locity dispersion, signal-to-noise ratio, and combination 
of stellar templates. To construct the simulated syn- 
thetic spectra we start with the best-fit linear combi- 
nation of stellar templates for three different regions in 
NGC 628 chosen to sample an old stellar population with 
a spectrum dominated by K giants (taken from the cen- 
tral part of the galaxy), a bright star forming HII region 
on a spiral arm, and a region showing a post-starburst 
spectrum dominated by A type stars. The three simu- 
lated spectra are convolved with Gaussian LOSVDs hav- 
ing o-* >tme =10, 25, 50, 100, 200, 300, and 500 km s _1 . 
Gaussian noise is artificially added to the spectra creat- 
ing 40 Monte Carlo realizations of each at S/N levels per 
resolution element (FWHM) of 25, 50, 100, and 200. 

We fit the synthetic spectra using the same methods 
used for the data and described in the previous section. 
Even though the simulated spectra do not contain emis- 
sion lines we mask the same regions which we censored 
when fitting the science data in order to match the avail- 
able spectral information. In an attempt to account for 
the effects of template mismatch we remove the tem- 
plates used to create the synthetic spectra from the list 
of templates used to fit them. The large number of tem- 
plates used (~ 50) allows us to do this without the worry 
of a particular stellar type not being represented in the 
remaining template subset. We present results for two 
fits, one using the full available spectral range, and one 
limited to the 5100-5300 A region around the Mg I b 
feature. 

The results of the simulation are presented in Figures 
llOl and fTTl In Figure [10] we compare the input and recov- 
ered values for the velocity dispersion at different S/N 
ratios. In Figure [TT] we present the error in the recovered 
radial velocity as a function of S/N for all the simulated 
spectra with cr* < 100km s _1 (typical of spiral galaxies). 
Results for fits using both the full spectral range and the 
limited 5100A-5300A window are reported. 

When using the full spectral range the recovered stel- 
lar velocity dispersions are overestimated by ~ 10% at 
f*. true > 100 km s^" 1 , become underestimated by ~ 30% 
below that value, and then become largely overestimated 
for values of cr* itrue which are about a factor of 10 smaller 
than the instrumental resolution. In the case of the full 
spectral range fit this trends are almost independent of 
the S/N ratio. As can be seen in the right panel of Fig- 
ure [TQ1 limiting the spectral window to the Mg I b region 
removes most of these systematic deviations at the cost 
of larger random error-bars on the recovered cr* values. 
Systematic deviations are still seen at low S/N and low 
c*,truc. Balmer lines play a significant role in introducing 
the systematic deviations seen in the full spectral range 
case as is confirmed by the fact that the effect is stronger 
for the "star-forming" and "post-starburst" spectrum, 
than it is for the "old stellar population" spectrum. No 
significant difference is seen on the recovered velocity dis- 
persions for the three different types of spectra once the 
wavelength range is limited to the Mg I b feature. 

An effect of limiting the spectral window is an increase 
in the random error in the recovered value of cr*. This 
increase is expected because of the high EW of the re- 
jected Balmer absorption lines. The error in cr* for each 
spaxel in the data-cube is computed by linear interpo- 
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Fig. 10. — Top panel: Comparison between the measured (o^/it) 
and the input (cr*,true) velocity dispersion from fits to Monte Carlo 
realizations of simulated spectra. Colors indicate the S/N of the 
synthetic spectra. Groups of datapoints for different values of 
a*,true are labeled in the upper left panel. Filled circles and error- 
bars indicate the mean and standard deviation of the recovered 
values from 40 Monte Carlo realizations. Points are shifted by 
small amounts along the horizontal axis for clarity. Bottom panel: 
Ratio between the recovered and input velocity dispersions. Dot- 
ted and dashed horizontal lines mark 20% (5<r measurement) and 
33.3% (3(r measurement) deviations respectively 



lation of the measured errors in the simulation to the 
spaxel's S/N and er„. 

The radial velocity is recovered with no systematic de- 
viations. Considering all simulated spectra with a n.^ruc < 
100 km s _1 (typical of the disks of spiral galaxies) and 
at all the S/N ratios probed, we measure mean offsets 
of —3 ± 11 and —2 ± 15 km s" 1 for the fits using the 
full spectral range and the 5100A-5300A window respec- 
tively. The error in the recovered radial velocity depends 
on S/N as shown in Figure [TT] and it is reduced by 
about a factor of two at the typical S/N ratio of the 
data by using the full spectral range including the strong 
Balmer lines. The improvement is even more dramatic 
at low signal-to- noise (S/N < 50), where the Mg triplet 
starts to get lost in the noise while the Balmer lines can 
still be easily seen in the spectra. Errors in for each 
spaxel are computed from the S/N by linear interpola- 
tion of the simulation results (black line in Figure Hip , 
and adding in quadrature the wavelength solution r.m.s 
of 6 km s" 1 (§4.2). 

5.3. Emission Line Fluxes and Ionized Gas Kinematics 

After measuring the LOSVD of each spaxel we use 
GANDALF to fit the spectrum including the emission 
lines. GANDALF recomputes the weights given to the 
different stellar templates at the same time of adding 
Gaussian profiles to model the contribution from emis- 
sion lines. We attempt to fit all the transitions listed in 
Table© 

During this fit the stellar LOSVD computed in the 
previous section is held fixed but the emission lines ve- 
locity, velocity dispersion, and amplitude, as well as the 
weights given to the stellar templates are free parame- 
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Fig. 11. — Error (standard deviation) in the recovered stellar 
velocity for fits to simulated spectra with o+^true < 100km s — 1 as a 
function of signal-to-noise ratio. Filled and open circles correspond 
to fits using the full spectral range and fits limited to the 5100A- 
5300A window respectively. 

ters. While in principle we could fit the kinematics of 
different lines independently, this becomes very hard for 
faint transitions detected at low S/N. Therefore we tie 
the kinematics of all emission lines to a common set of 
parameters (v gas , o~ gas ) during the fit. This ensures that 
the kinematic parameters are dominantly constrained by 
the brightest emission lines in the spectrum (typically 
Ha, [OIII]A5007, and [OII]A3727). In this second iter- 
ation we also use additive and multiplicative Legcndrc 
polynomials to match the continuum shape. 

Figure [TJ] presents the observed and best-fit spectrum 
of a randomly selected spaxel in the data-cube. This 
spectrum provides a good representation of the quality 
of our fits. The observed spectrum is shown in blue, with 
photometric errors marked by the cyan envelope, and to- 
tal errors (photometric plus systematic) shown as a green 
envelope. The solid red line is the best-fit stellar plus 
emission line spectrum and the dotted red line shows the 
best-fit stellar spectrum alone. The vertical cyan bands 
represent regions masked around sky line residuals and 
the ends of the blue and red setup spectra. Also shown 
are zoomed in windows around relevant spectral features. 
Overall GANDALF produces fits to the VENGA spectra 
of excellent quality. 

5.4. Errors in Emission Line Parameters and 
Kinematics 

GANDALF returns emission line parameter errors 
based on the covariance matrix computed during the 
non-linear Levenberg-Marquardt fitting. Therefore, the 
quality of this error estimates depends critically on the 
uncertainty in the observed spectrum being properly esti- 
mated. In particular for low EW emission lines detected 
in very high S/N spectra, for which the formal photomet- 
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Fig. 12. — Observed and best-fit spectrum of a randomly selected spaxel in the NGC 628 data-cube. The spectrum has a S/N ratio of 
87 in the continuum per spectral resolution element. The observed spectrum (solid blue) is shown with photometric errors (cyan envelope) 
and total errors (photometric plus systematic, green envelope). Also shown is the best-fit stellar plus emission line spectrum (solid red), 
and the best-fit stellar spectrum alone (dotted red). Vertical cyan bands represent regions masked around sky line residuals and the ends 
of the blue and red setup spectra. Also shown are zoomed in relevant spectral features. 
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Fig. 13. — Maps of NGC 628 in different emission lines. From 
top to bottom: Ha, H/3, [OII]A3727, [OIII]A5007, [NII]A6584, and 
[SII]A6716. Black contours show lines of constant galactocentric 
radii spaced by 1 kpc, with the exception of the innermost contour 
at a radius of 0.5 kpc. Spaxels contaminated by foreground MW 
stars (see §6.1) are masked. 

ric errors are small, the inclusion of the systematic uncer- 
tainties associated with template mismatch (described in 
§5.1) is crucial in order to obtain adequate emission line 
flux errors. If this floor of systematic uncertainty is not 
included, errors in the measured emission line parameters 
can be significantly underestimated. 

Figure Q~3] shows maps of the flux of the Ha, H/3, 
[OII]A3727, [OIII]A5007, [NII]A6584, and [SII]A6716 lines 
across NGC 628. We only show spaxels in which the line 



is detected with a significance higher than 3(7. The Ha: 
line is detected at 5(7 over the full data-cube. The neb- 
ular emission clearly traces the two main spiral arms of 
the galaxy where ongoing star formation gives rise to 
prominent HII regions. In the inter-arm regions we still 
detect significant amounts of Ha emission although with 
a surface brightness that is one to two orders of mag- 
nitude fainter than in the arms. A large fraction of this 
emission arises from the diffuse ionized gas component of 
the galaxy's ISM and is enha nced in low ionization ratios 
like [SHI /Ha and [Nil] /Ha (|Mathidl200CH IHaffner etHI 
[200a iBlaric et all 12009). Also clearly seen in Figure [H 
is a circum-nuclear star forming ring with a diameter of 
~ 0.5 kpc and a Ha hole in the center of the galaxy. 
These features and their physical origin will be discussed 
in a future paper. 

For each transition in Table |H] we report the median 
S/N over all spaxels in the data-cube and the fraction of 
the observed area in which the emission line is detected 
at 5<7 and 3<7. Our data on NGC 628 reaches a 5<r line 
flux limit per spaxel (4 arcsec 2 ) of ~ 2 x 10~ 17 erg s _1 

_2 

cm . 

As mentioned above, we detect Ha at 5cr over the 
full data-cube with a median S/N — 47. Other 
transitions are usually detected at lower significance 
than Ha. Individual members of the [NII]AA6548,6583 
and [SII]AA6717,6731 doublets, as well as H 7 , H/3, 
[OIII]A5007, and the blended [OII]AA3726,3729 doublet 
are detected at 3a over more than ^90% of the observed 
area, and the [OIIIJA4959 line is detected at 3a over 58% 
of the area. 

5.5. Nebular Extinction from Balmer Decrement 

For the typical extinction regimes found in spiral galax- 
ies the Ha to H/3 flux ratio (a.k.a Balmer decrement) 
provides a good estimate of the amount of dust extinc- 
tion in the region s giving rise to th e nebular emission 
(see discussion in I Blanc et al.l [2009) . Assuming an in- 
trinsic Ha/H/3 ratio of 2.87 for Case B recombination in 
n = 10 2 cm" 3 , T = 10 4 K gas (|Osterbrock fc Ferlandl 
2006), the observed ratio provides the amount of dust 
reddening by means of the following relation 

-2.5 log ( ra^) 

E(B -V)= , V - f J (4) 

k(Ah q ) - fc(AH/3) 

where [Ha/H/3] f, s is the observed line ratio and k(X) is 
the extinction law. We assu me a foregr ound MW extinc- 
tion law as parameterized bv lPel (I1992D. S MC and LMC 
laws were tested (also using the lPeilll992l parametriza- 
tion), and no significant change was observed in the de- 
rived extinction values (these 3 extinction laws are prac- 
tically identical at optical wavelengths). 

Figure [T4l presents a map of E(B — V) across the disk 
of NGC 628. Errors in the measured values of E(B — V) 
are computed by propagating the uncertainty in the Ha 
and H/3 line fluxes. In Figure [14] we have masked re- 
gions in which the error in E{B — V) is larger than 0.2 
mag, which roughly amounts to an uncertainty of ~ 50% 
in the extinction correction factor at the wavelength of 
Ha. Regions of high extinction follow the inner rims 
of the spiral structure (as traced by Ha), and a visua l 
comparison to CO maps of NGC 628 (|Lerov et al.ll2008l) 
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shows almost perfect coincidence between high dust ex- 
tinction and high molecular gas surface density regions. 
Also presented in Figure [14] is a map of Ha after correct- 
ing for dust extinction. The Balmer decrement method 
implies a median extinction across all regions with good 
E(B — V) measurements of Ay = 1.06 mag (assuming 
R v = A V /E{B - V) = 3.1). 

In Figure [T5] we compare the VENGA Balmer decre- 
ment dust extinction measurements to a compilation of 
HII region measurements in the literature. The VENGA 
datacube of NGC 628 covers 2 HII regions from the 
catalog of iRosales-Ortega et al.l ([201 If ) (built from the 
PINGS IFU datacube of the galaxy ) 4 HI I regions from 
the multi-slit stud y of |Gusev et all (I2012D . and two HII 
regions from the iMcCall et al.l (I1985D sample, one of 
which was also studi ed by Ferguson et al.l (fl998). The 
two regions from the IMcCall et al.1 (|1985[ ) sample were 
also detected by PINGS so the total number of unique 
comparison HII regions in 24. In our literature search 
we only considered regions with statistically significant 
extinction measurements in the respective studies. 

Beacuse of the large spatial resolution of our data the 
comparison to seeing limited slit spectroscopy measure- 
ments is challenging. Even the PINGS data has a spatial 
resolution which is a factor two smaller than ours. For 
the VENGA measurements we use the average extinc- 
tion in a 5.6" diameter aperture centered in each region, 
expecting the nebular emission in that spatial resolution 
element to be dominated by the HII region. Consider- 
ing the differences in methodology the comparison shows 
a reasonable agreement between the values measured in 
VENGA and those measured by previous authors. The 
observed scatter is consistent with the scatter seen in ex- 
tinction among different st udies. For the one HI I regio n 
which is common to the IRosales-Ortega et al.l (| 2 1 If ) . 
IMcCall et all (|1985| ). and iFerguson et all (|1998T ) studies 
the authors measure Ay values of 0.38, 0.56, and 1.66 
mag respectively while the VENGA data yields an inter- 
mediate value of 1.27. 

6. FITTING THE STELLAR AND IONIZED GAS VELOCITY 

FIELDS 

In this section we present the methods used to fit the 
stellar and ionized gas velocity fields of NGC 628. We are 
interested in measuring the galaxy's rotation curve using 
both the stellar and gaseous component with the goal 
of understanding any differences observed between both 
velocity fields. We propose a new method to measure 
the inclination of face-on disk galaxies which is based on 
matching the gas and stellar rotation curves by means of 
the asymmetric drift correction. 

6.1. Removal of Foreground Stars 

In Figure [TH] we present the stellar and ionized gas 
velocity field of NGC 628 resulting from the fits described 
in §5. Foreground Milky Way stars can be clearly seen in 
the stellar velocity field as sources of low velocity along 
the line of sight. This is evident by looking at the SDSS 
r-band isophotes overploted on the top panel of Figure 
[Tin Guided by the SDSS broad-band image, we flag these 
low velocity spaxels in the data-cube as spectra which 
are significantly contaminated by foreground stars. We 
mask all regions within 5" of these flagged pixels. The 
resulting stellar and ionized gas velocity fields after the 
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Fig. 14. — Upper Panel: Map of the dust reddening E(B — V) 
computed using Equation 4. Regions where the uncertainty in 
E(B — V) is larger than 0.2 mag have been masked. These masked 
regions correspond mainly to low surface brightness parts of the 
galaxy where the Ha to H/3 ratio cannot be measured with high 
S/N. Lower Panel: Map of the Ha emission line flux corrected for 
the effects of dust extinc tion using a MW extinction law. Black 
contours are as in Figure [13] Spaxels contaminated by foreground 
MW stars (see §6.1) are masked in both panels. 
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Fig. 15. — Comparison between the dust extinction (Ay) derived 
from the Balmer decrement in VENGA and measru ements found 
in the literature for a sample of HII regions in the IMcCall et al.l 
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1T985I); Ferguson eFaD GMD; fRosalcs-O rtega et al.l 1 1201 U) and 
I Gusev e t al. (2012j) samples. Dashed lines indicate a 50% devi- 
ation in the dust correction. 



masking of foreground stars are shown in the second and 
bottom panels of Figure [TB] Even at the low inclination 
of NGC 628 (see below for a discussion on the actual 
value) the VENGA data allows us to clearly detect the 
rotation of both the stellar and gaseous components of 
the galaxy. 

6.2. Kinematic Position Angle and Systemic Velocity 

We use the s ymmetrization m e thod described in 
Appendix C of iKrainovic et al.l (|2006[ ) and imple- 
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Fig. 16. — First Panel: Map of the stellar LOS velocity (ti*) in 
NGC 628. Black contours mark SDSS r-band isophotes as in Figure 
[9] Gray dots mark spaxels flagged as contaminated by foreground 
stars. Second Panel: Clean map of the stellar LOS velocity where 
both foreground stars and regions with an error in larger than 
15 km s have been masked. Third Panel: Map of the LOS stellar 
velocity dispersion (cr* = oxos) meas u re d from the fit to the Mg I 
b region. Lower Panel: Map of the ionized gas LOS velocity (v gas ) 
as measured from the simultaneous fit to all emission lines listed 
in Table [6] Regions having an error in v gas larger than 15 km s 
have been masked. The zero-velocity line for the measured 21° 
PA is marked by the diagonal solid black line in the bottom three 
panels. 

mented by Michele Cappellari in the IDL routine 
FIT_KINEMATIC_PAE!I to measure the kinematic po- 
sition angle (PA) and systemic velocity (v sys ) from both 
maps. In brief, for a given set of values (PA, v sys ), the 
method assumes mirror symmetry in the velocity field 
and a symmetrized velocity field is constructed by av- 
eraging the absolute values of the velocity at mirroring 
points on opposite sides of the perpendicular to the line of 
nodes. The best-fit PA and v sys are those which minimize 
the residuals between the observed and symmetrized ve- 
locity fields. 

For the kinematic PA we obtain consistent values 
of 20 ± 5° and 22 ± 6° for the stellar and ionized 
gas velocity fields respectively. This value is in agree- 
ment with the RC3 J3-ba nd photometric PA of 25° 
(jde Vaucouleurs et al.lll99lT ) , but is inconsi stent with the 
2MA SS if-band reported value of 87.5° (jJarrett et al.1 
2003). The kinematic PA of the ionized gas velocity 



field in the inner parts of the disk of N GC 628 has 
been m easur ed independ e ntly in the past bv lDaigle et al.l 
(|2006T) and iFathi et al.l (|2007D using Fabry-Perot Ha 
spectroscopy. These two studies measured inconsistent 
values of 26.4 ± 2.4° and 15 ± 5°. The value derived from 
the VENGA data sits right in between these two mea- 
surements and we consider it to be more reliable as it 
is based on deeper data, using the kinematics of multi- 
ple emission lines, and shows good agreement with the 
k inematic PA of the stellar c omponent . 

iKamphuis fc Briggsl (|1992P ) used VLA HI 21 cm obser- 
vations to reveal the presence of an extended HI disk 
in NGC 628 reaching out to a galactocentric radius of 
^30 kpc and showing complex morphology and kine- 
matic structure. The HI velocity field shows a sharp 
transition in both PA and inclination around a galac- 
tocentric radius of 7' (~ 18 kpc) believed to be caused 
by the presence of a warp in the outer HI disk. The 
VENGA IFU data does not sample the disk out to such 
large galactocentric radii but the stellar and ionized gas 
kinematics are consistent with the HI kinematics within 
the inner disk of the galaxy. The HI velocity field has 
a kinematic PA of ~ 25° wit hin the central 4' (Figure 6 
in IKamphuis fc Bri ggs (1992)), consistent with our mea- 
sured value. 

The symmetrization procedure also yields the sys- 
temic velocity of the galaxy. For the gaseous compo- 
nent we obtain a value of 650 ± 10 km s , in agree- 
ment wjHi__ft£jHn_2l£m derived values of 657 ± 1 km 
s" 1 (IKamphuis fe Briggsl 119921 ) and 659 ±3 km s" 1 
(I Walter et al.ll2008ri . but in stark disagreement with the 
value of 683 ± 3 km s" 1 derived by IFathi et al.1 (|2007P ) 
from the Ha velocity field. This inconsistency proba- 
bly reflects a syste matic error i n the p arameters of the 
kinematic model in IFathi et al.l (|2007f ) who also report 
a decreasing rotation curve within the central kiloparsec 
of the galaxy, feat ure which is not observed neither by 
iDaigle et al.l(|2006l ) nor in this study. The VENGA stel- 
lar velocity field yields a consistent systemic velocity of 
651±10 km s _1 also agreement with the HI 21cm derived 
values mentioned above. 

6.3. Harmonic Decomposition Modeling 

In order to measure the rotation curve of the galaxy we 
model the observed stellar and ionized gas v elocity fields 
using the harmonic decomposition method (iFranx et ail 
IT99llWong et al.ll2004l:rKrainovic et al.ll2006D . Inpartic- 
ular we use the IDL KINEMETR M3 package to perform 
the fitting ()Krainovic et al.ll2006|) . We fit the azimuthal 
dependance of the LOS velocity along a set of ellipti- 
cal annuli at different radii along the semi-major axis 
of NGC 628. The annuli semi-major axes are spaced fol- 
lowing the combinat ion of linear and l ogarit hmic increase 
factors described in iKrajnovic et al.l (|2006[ ) with a min- 
imum radius of 6" . For each ellipse the LOS velocity is 
fit as 

3 

v(R,ip) = v sys + ^[A n (i?)sin(m/>) + B n (R)cos(mp)] 

71=1 

(5) 



17 http:/ /www-astro. physics. ox. ac.uk/~mxc/idl/ 
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where v sys is the systemic velocity measured indepen- 
dently for each component in the previous section, R 
is the ellipse semi-major axis, and ip is the eccentric 
anomaly angle. The outermost radius at which the fit 
is performed is set by the requirement that data is avail- 
able for 70% of the points along the ellipse. For NGC 
628, this translates in an outermost radius of 60" (~ 2.5 
kpc). _ 

During the fitting, the geometric parameters of the el- 
lipses (center, PA, and inclination) are held fixed. We 
use the average of the two globally derived kinematic 
PAs for stars and gas of 21° (§6.2), and the central co- 
ordinates from NEEpT. As discussed below, the choice of 
an inclination value for NGC 628 is far from trivial. 

Using photometric axis ratios (q — b/a) to com- 
pute t he inclination of NGC 628 by adopting the iTullvl 
(|1988T> method we derive i = 28° (q = 0.91) for the 
RC3 B-band measurements ()de Vaucouleurs et al.lll99lT) 
and i = 34° (q = 0.86) for the 2MASS if-band pho- 
tometry (| Jarrett et al.l I2003D . The detailed photomet- 
ric decomposition of iFisher k Drorvl (|2008D implies an 
inclination of i = 33° (q = 0.87) for the disk com- 
ponent. Given the projected circular velocities mea- 
sured at l arge radii in both the 21cm HI and Ha veloc- 
ity fields | Kamphuis_&BriggsJ Il992t IDaigle et all 120061: 
IFathi et all 120071) . such inclinations would imply a de- 
projected velocity of Vfi a t ~ 50 km s _1 for the fiat part 
of the rotation curve. Such a small circular velocity 
is wildly inconsistent with expectations fr om the Tully- 
Fisher (TF) relation (jTullv k Fisherlll977l) for NGC 62 8. 

Using the if -band TF relation from lVerheiienl (<2001[F 1 
and the 2MASS if-band absolute magnitude from Ta- 
ble [T] we obtain u/; a t,TF = 149 ± 3 ± 8 km s _1 , where 
the random error comes from the uncertainty in the dis- 
tance and apparent magnitude of the galaxy and the 
systematic uncertainty comes from the 0.26 mag scat- 
ter in the adopted TF relation. The inconsistency be- 
tween this value and that derived from the photomet- 
rically measured inclination is not surprising given that 
the uncertainty in i explodes as one approaches a nearly 
face-on configuration. This is also the case for kinemat- 
ically derived inclinations. In fact, letting the inclina- 
tion be a free parameter for the ellipses in our harmonic 
decomp osition fits yields an unrealistic median value of 
i = 45°. Daigl e et al.l (|2006) derives a kinematic inclina- 
tion of 21.5 ± 4.5° from their tilted ring analysis of the 
Ha velocity field. For an excellent discussion on the un- 
certainties assoc iated with measuring th e inclination of 
disk galaxies see IBershadv et al.l (|2010b|) . 

Given the unreliability of both photometric and kine- 
matic inclinations for face-on galaxies, most studies of 

NGC 628 have assumed a value of i = 6.5° (e.g. 

I Kamphuis k Briggslll992t IDaigle et al.ll2006HF athi et al.l 
l2007HHerbert-Fort et allboiOl) . This value was adopted 
by IKamphuis k Br iggs (1992|) so the flat part of the HI 
rotation curve would reach an arbitrary chosen value of 
Vfiat = 200 km s _1 . A notable exception is the work 



19 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 

20 We use the parameters from th e fit to the i r "RC /FD Sample, 
Excluding NGC 3992" in Table 4 of I Verheiier] [120011) 



of iHerrmann et all (|2008l) who fit an exponential model 
to the THINGS 21 cm HI velocity field of NGC 628 to 
obtain an inclination of i — 9°. 

Under the assumption that NGC 628 exactly follows 
the TF relation, a better value for the inclination would 
be %tf — 8.7 ± 0.7° (assuming an error of 10 km 
s -1 for the value o f v fiat measured at large radii by 
IKamphuis k Briggsl (|1992l ) and propagating the errors 
in the apparent magnitude and distance of NGC 628 
as well as the scatter in the TF relation). This so 
called "Tully-Fisher inc lination" method is discussed in 
IBershadv et al.l (|2010bD and provides more reliable esti- 
mates of the inclination for face-on systems than photo- 
metric or kinematic techniques do. 

In the following section we propose a new alternative 
method to derive the inclination of face-on disk galax- 
ies by demanding agreement between the ionized gas 
and stellar rotation curves after applying the asymmetric 
drift correction to the stellar component. 

6.4. Rotation Curve and Asymmetric Drift Correction: 
A New Method to Measure Inclinations 

Because of the collisional nature of gas, as opposed 
to the collisionless nature of stars, the cold ISM of disk 
galaxies is expected to be confined to a relatively thin 
disk supported mainly by rotation. Although pressure 
support in the ISM can become important under cer- 
tain circumstances (e.g. high turbulence, strong pres- 
sure gradients, etc.), the small and relatively constant 
Ha velocity disper sion of a Ha ~ 17 km s" 1 measured by 
IFathi et al.l (|2007t ) across the disk of NGC 628 implies 
that the gaseous disk of this galaxy is almost fully sup- 
ported by circular motions, with random motions only 
having a minor effect. The importance of pressure sup- 
port against gravi ty in gaseous galac t ic dis ks is exten- 
sively discussed in lDalcanton k Stilpl (120101) . Using the 
flat dispersion pro file of Fathi et al. ( 20071) and applying 
equation 11 from IBalcanto n k Stilpl ( 20101 ) we compute 
a correction to the ionized gas rotation curve of NGC 628 
of ~ 2 km s 1 . This is many times smaller than the la 
uncertainty in the rotation curve itself. Given the small 
magnitude of this effect we decide not to apply these cor- 
rections and neglect the effect of pressure support in the 
gaseous disk of NGC 628. Therefore we consider the ion- 
ized gas rotation curve to roughly follow the dynamics 
of a test particle under the gravitational potential in the 
mid-plane of the galaxy. 

Stellar disks can also experience support against grav- 
ity provided by random motions. As in the case of 
pressure support in gaseous disks, stellar random mo- 
tion support translates in a measured circular velocity 
which is smaller than the actual circular velocity at which 
a test particle would orbit in the gravitational poten- 
tial of the galaxy. This effect, first noted as a differ- 
ence in orbital velocity for stellar populations of differ- 
ent ages (and therefore different velocity dispersions) in 
the solar neighborhood is kno wn as "asymmetric drift" 
(jStromberell922D . As noted bv lDalcanton k Stilpl (|2010f > 
the equations describing it are similar to those describ- 
ing pressure support in gaseous disks but the physical 
phenomena are fundamentally differ ent. 

Star ting from the Jeans equation s (|Jeansll 19151 ) and fol- 
lowing iBinnev k Tremaind (|1987l) , for a stellar axisym- 
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metric disk in steady state, under the assumption that 
the "epicyclic approximation" holds and that the stellar 
velocity ellipsoid (SVE) is aligned with the disk in cylin- 
drical coordinates and its shape is independent of the 
vertical distance to the plane, the quadratic difference 
between the true and the observed circular velocities is 
given by 



a% f dln(u ) 
2 \dln(R) 



1 



J d\n(R) 



(6) 



where v c is the circular velocity of a test particle in the 
potential, is the observed circular velocity, an is the 
velocity dispersion in the radial direction, and v is the 
number density of stars. While and its logarithmic 
derivative with radius can be directly measured from the 
harmonic decomposition fits described in §6.3, calculat- 
ing the last term in the perenthesis in Equation 6 requires 
making a series of assumptions which we describe below. 

In order to apply the asymmetric drift correction we 
assume a constant scale height and V-band mass-to-light 
ratio for the disk of NGC 628. These assumptions imply 
that v(R) is proportional to the stellar mass surface den- 
sity and that the latter is well traced by the surface 
brightness profile of the galaxy. In particular we use the 
best-fit disk plu s pseu do-bulge V-band decomposition of 
iFisher fc Drorvl (|2008D to trace the radial dependance of 
v{R) when calculating the last term in Equation 6. Note 
that knowledge of the V-band mass-to-light ratio is not 
necessary as it cancels out in the logarithmic derivative. 

Using a measurement of the stellar velocity dispersion 
along the line of sight (olos(-R)) we can estimate an(R), 
provided our assumptions regarding the shape and orien- 
tation of the SVE. The third panel in Figure [TBI presents 
the map of the stellar velocity dispersion (er* = ctlos) 
measured from the VENGA data in §5.1. In the central 
regions of NGC 628 we see cr» falling from a central value 
of ~ 60 km s" 1 . Figure [T71 presents the radial profile of 
clos- To measure the radial profile we only considered 
spaxels with high enough S/N as to provide a reliable fit 
of the LOSVD (i.e. those which are not masked in Figure 
I16[) . Red circles mark the median ctlos m 0.5 kpc wide 
radial bins. Error bars show both the standard deviation 
in each radial bin (red) and the median error in for 
each spaxel (black). 

Also plotted in Figure [T7] is the central stellar ve- 
locity dispersion value of a*(R = 0) = 72.2 ± 7.8 km 
s- 1 from HyperLedcE] (blue circle). The HYPERLEDA 
value is an average over six literature mea surements and 
is dr iven by an outlier at 116 km s^ 1 (di Ncll a et al.l 
I1995D - Rejecting this outlier the HYPERLEDA aver- 
age comes down to 64 km s _1 , in even better agree- 
ment with our measurements (blue triangle in Figure 
I17[) . The velocity dispersions we measure in the cen- 
tral regions of the galaxy are also consistent with the 
nuclea r value of 58.0 ± 8.6 km s _1 measured bv lHo et"al] 
( 2009) (green circle) . Our measurements also agree with 
the central d i spersi on of 54 ± 10 km s _1 measured by 
iGanda et all (|2006ft using SAURON IFU spectroscopy. 
The radial cr„ profile measured in this last study is shown 
as a dashed black line in Figure [T7] While our mea- 

21 http://leda.univ-lyonl.fr/ 
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Fig. 17. — Radial profile of a* = <7los- Red circles mark the me- 
dian cr* for all spaxels in 0.5 kpc wide radial bins. Error bars show 
both the standard deviation in each radial bin (red) and the me- 
dian error in <r» for each spaxel (black) . The cent r al ste llar velocity 
dispersion values from HyperLeda and Ho et al. (2009) arc shown 
in blue and green r espectively. The dispersion profile measured by 
IGanda ct al. (2006) using the SAURON IFU spectroscopy is shown 
as a black dashed curv e. The velocity d i spers ion profile measured 
from PN kinematics by Herrmann ct al. (2008) is also shown (light 
blue circles, dashed line). The solid black curve shows the best-fit 
dispersion profile to the VENGA data used to calculate the asym- 
metric drift correction. 



sured values in the center of the g alaxy are consistent 
with the ones in IGanda et al.l (|2006l ) we find no evidence 
to support their claim of an increasing velocity disper- 
sion with radius in NGC 628. At large r radii Figure [TTI 
shows the dispersion profile measured bv lHerrmann et all 
(2008) using the kinematics of a large sample of plan- 
etary nebulae (PN). Th e dashed line sh o ws th e best 
fit dispersion profile from iHerrmann et al.1 ()2008[ ) which 
follows an exponential with a scalength aproximatcly 
equal to two times the photometric scalelength of the 
galaxy. The latter agrees with theoretical expecta- 
tions for constant scaleheight disks in which a z (R) oc 
^Tidisk(R), with the proportionality constant depending 
on the detailed vertical profile of the disk (exponential, 
sech or sech , Ivan der Kruidll988t iBershadv et alJl2010t 
Ivan der Kruit fc Freeman! 120 111 and references therein). 

The 145 km s _1 spectral resolution at Mg I b in the 
VENGA VIRUS-P data makes the measurement such 
low values of the velocity dispersions very challenging. 
The large error bars in Figure [T7] indicate that for single 
spectra individual measurements of low dispersion val- 
ues (~ 50 km s _1 ) are extremely noisy even at the high 
S/N ~ 100 of our data. Nevertheless, by averaging over 
many spaxels in the data-cube <Jlos we seem to recover 
the velocity dispersion profile of the galaxy at least out 
to R ~ 4 kpc where the VENGA measurements flatten 
around a value of ~ 40 km s _1 . Comparison to the dis- 
persion profile derived from PN kinematics shows that 
this flattening is not real and most likely corresponds 
to a systematic floor caused by uncertainties associated 
with template missmatch and deviations from the as- 
sumed Gaussianity in the instrumental line spread func- 
tion most likely set this limit. This is not surprising 
given that it happens at about a fourth of the instru- 
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mental resolution and that a similar behavior is seen in 
the simulations presented in §5.2. 

We fit an exponential dispersion profile to the VENGA 
data in the 1.5 < R < 4 kpc region in order to avoid 
regions in which the stellar dynamics are contaminated 
by the central spheroi dal component (see discussion in 
iBershadv et al.l l2010bD . We fix the scale length of the 
profile to t wo times the scale leng th of the disk compo- 
nent in the iFisher fc Drorvl (|2008f ) decomposition. The 
best fit profile is shown as a solid line in Figure [TTl im- 
plies a central velocity dispersion of ( 7t ; or(0) = 62 km s 
(consistent w i th the HYPERLEDA, |Hg_et_al] (|2009fl and 
iGanda et al.l (|2006[) value s ) and is in excellent agreement 
with the lHerrmann et al.l (|2008T ) best fit profile and with 
theoretical expectations for the dynamics of stellar disks. 
This itlos(^) profile is used to derive aB.(R) in order to 
calculate the assymetric drift correction using Equation 
6. 

Follo wing the method outlined in Bcrshadv et al.l 
(2010b) we assume the following shape for the SVE in 
order to deproject the LOS velocity dispersion profile 
and compute <jr(R): 
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where the first equatio n comes from expectations from 
the solar neighborhood (|Binnev fc Merrifi cld 199 8}) and 
observations of nearby galaxies ()Bershadv et al.l 12010b) 
and the second comes from adopting the epicyclic ap- 
proximation. This allows us to compute all the terms 
in Equation 6 and apply the asymmetric drift correction 
to the stellar velocity field of the galaxy. We use Monte 
Carlo simulations to calculate the error in the asymet- 
ric drift correction. In the error calculation we include 
the uncertainty in the observed rotation curve which is 
propagated to compute the uncertainty in the logarith- 
mic radial derivative of the rotation curve and the value 
of j3, an error of Aa = 0.15 ( jBershady et al.lfeOlObl ). and 
an error in ctlos(O) of 10 km s _1 . 

In the upper panel of FigureHUwe present the VENGA 
rotation curve of NGC 628 for both stars (black circles) 
and ionized gas (red circles), as measured from the first 
order cosine term of the harmonic decomposition model 
Vcirc(R) = Bi(R)/sin(i). We show measurements of the 
rotation curve for two different assumed inclination val- 
ues of i = 28° (derive d from the RC3 axis ratio s, open 
circles) and i — 6.5° (jKamphuis fc Briggslll992l , filled 
circles). The horizontal dashed line marks the expecta- 
tion from the TF relation Vfi a t,TF = 149 km s _1 (§6.3). 
The high inclination implied by the photometric axis ra- 
tio dramatically underestimates the amount of rotation 
expected in NGC 628. This highlights the lack of reli- 
ability of geometrically measured inclinations in face-on 
systems. 

Independent of the assumed inclination, the observed 
stellar circular velocity (before applying the asymmetric 
drift correction) is systematically lower than the observed 
circular velocity of the gaseous component. This implies 
the presence of significant kinematic support arising from 
random motions in the stellar disk of NGC 628. For the 
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Fig. 18. — Top Panel: Stellar (black) and ionized gas (red) ro- 
tation curves for NGC 628 under the assumption of two different 
inclination values of 6.5° (filled circles) and 27.8° (open circles). 
Open squares show the stellar rotation curve after applying the 
asymmetric drift correction for the i = 6.5° case. The horizontal 
dashed line marks the asymptotic rotation velocity of 149 km s 
expected from the Tully-Fisher relation. The size of the VENGA 
PSF FWHM is marked by the dotted vertical line. Lower Panel: 
Same as above, but for the best-fit asymmetric drift inclination 
(iAD = 8.7° ± 0.7). The inset shows the \ 2 distribution for the 
range of inclination values explored during the fit. The vertical 
dotted line in the inset marks the best-fit inclination. 



i = 6.5° case the corrected stellar rotation curve after ac- 
counting for asymmetric drift using equation 6 is shown 
by the black open squares in the top panel of Figure 
[TBI Errorbars in the corrected stellar rotation curve are 
calculated using Monte Carlo simulations and including 
the uncertainty in all the quantities entering the correc- 
tion. For an assumed value of i = 6.5° the magnitude 
of the asymmetric drift correction is not sufficient to ac- 
count for the discrepancy between the observed stellar 
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and gaseous rotation curves. 

The above problem can be solved by changing the 
assumed value for the the inclination. We remind 
the reader that the i = 6.5° value typically used in 
the literature for NGC 62 8 was arbitrarily adopted by 
(jKamphuis fc B riggs 1991) to make the flat part of the 
rotation curve reach a velocity of 200 km s _1 . From 
the TF relation we have seen in §6.3 that the expected 
asymptotic velocity in the flat part of the rotation curve 
of NGC 628 is smaller than what was assumed by 
iKamphuis fc Briggsl (|1992l) with a value closer to 150 km 
s _1 . This asymptotic velocity implies an inclination of 
%TF = 8.7°. 

The lower panel in Figure [TSl shows the result of 
matching the gas and asymmetric drift corrected stel- 
lar rotation curves to obtain the inclination iad of the 
galaxy. The method consists of sampling a range of incli- 
nation values in order to find the inclination which min- 
imizes the residuals between both rotation curves. For 
each inclination value we measure the rotation curves 
by performing harmonic decomposition fits to the ve- 
locity fields, we deproject the LOS velocity dispersion 
profile to compute the asymmetric drift correction, and 
we compute a \ 2 statistic in order to quantify deviations 
between the gas and corrected stellar rotation curves. 

The inset in the lower panel of Figure [18] shows the 
\ 2 distribution as a function of inclination. The mini- 
mum x 2 corresponds to a best-fit value for the asymetric 
drift inclination of «ad = 8.7 ± 0.7° in perfect agree- 
ment with the %tf = 8.7 ± 0.7° value derived from the 
TF relation. The exact agreement in both the value and 
the error-bars implies that both methods produce con- 
sistent results with similar precisions, at least in the case 
of NGC 628. The TF and asymmetric drift methoda are 
highly complementary. The TF method requires knowl- 
edge of the observed rotation velocity in the flat part 
of the rotation curve at large radii, while the asymmet- 
ric drift method can be used even when only kinematic 
measurements are only available in the central regions 
of galaxies. Furthermore, the measurement of an inde- 
pendent asymetric dirft inclinations can allow for highly 
inclined systems to be included when measuring the TF 
relation. 

We propose the matching of the gaseous and asymmet- 
ric drift corrected stellar rotation curves as a new method 
to estimate the inclination of nearly face-on disk galax- 
ies. The fact that a very small difference of 2.2° in the 
adopted inclination can produce a catastrophic inconsis- 
tency between the observed gas and stellar kinematics 
as in the i = 6.5° case highlights the sensitivity of this 
method as tool to measure inclinations. 

7. SUMMARY AND CONCLUSIONS 

In this work we have presented the survey design, 
sample, observing strategy, data reduction, and spectral 
analysis methods for VENGA. Wide field integral field 
spectroscopy proves to be a powerful tool to study a large 
set of physical phenomena occurring in nearby galaxies. 
We characterized our sample of 30 nearby disk galaxies in 
terms of their stellar masses and SFRs, and showed that 
VENGA is a representative sample of massive (> 10 10 
M©) spiral galaxies in the local universe. A large range of 
morphological parameters (bulge-to-disk ratio, bar pres- 
ence, bulge Sersic index, inclination) is represented in the 



sample. The distance distribution of our targets implies 
a median spatial resolution in physical units of 383 pc. 
Our targets typically span a few arcminutes in angular 
size, and we tile the VIRUS-P IFU in order to provide 
a coverage out to 0.7 R25 for most galaxies. This im- 
plies that the VENGA data-cubes are among the largest 
optical IFU data-cubes ever constructed. In particular, 
the NGC 628 data-cube presented in this paper spans 
5.2' x 1.7'. 

The wealth of information produced by integral 
field spectrographs stresses the need for optimized and 
pipelined software tools to process and analyze the data. 
This becomes essential when this observational technique 
is used to conduct large surveys like the one presented 
here. We have presented the reduction and calibration 
pipeline used for the VENGA data. We have also de- 
scribed our spectral analysis pipeline which we use to 
extract stellar and gas kinematics as well as emission 
line fluxes, and we have discussed the methods used to 
estimate the uncertainties in all these parameters. When 
possible, we have adapted existing publicly available soft- 
ware to be used on the VIRUS-P data. We assessed the 
quality of the VENGA data obtained on NGC 628 and 
we find it to be excellent. VIRUS-P provides high S/N 
spectra for single fibers out to large galactocentric radii. 
Thanks to the ability of reconstructing broad-band im- 
ages from the IFU data and cross-correlating them with 
archival broad-band images of the galaxies in our sample 
we can achieve good astrometric precision (~ 0.1") and 
better than 10% flux calibration accuracy across the full 
wavelength range. 

Finally we have used an harmonic decomposition tech- 
nique to fit the stellar and ionized gas velocity fields 
in NGC 628. The observed kinematics imply that the 
gaseous disk of the galaxy (as traced by nebular emission 
from ionized gas) is almost fully supported by rotation 
while the stellar disk shows significant levels of dynamic 
support arising from random motions. This is evident by 
a ~ 50 km s^ 1 offset between the ionized gas and stellar 
rotation curves. 

We have shown that the above discrepancy in the rota- 
tion curves can be explained if one takes into account the 
effect of stellar random motions by applying an asymmet- 
ric drift correction. The agreement between the ionized 
gas rotation curve and the asymmetric drift corrected 
stellar rotation curve depends strongly on the assumed 
value for the inclination. We take advantage of this 
dependance and propose this as a method to estimate 
the actual inclination of nearly face-on disk galaxies for 
which classically used geometrically derived inclinations 
are extremely unreliable. Applying this method we find 
a revised asymmetric drift inclination for NGC 628 of 
iad = 8.7 ± 0.7°, in excellent agreement with the ex- 
pected value of i — 8.7±0.7°, which comes from assuming 
that NGC 628 satisfies the local TF relation. 

The main purpose of this work is to present the 
VENGA survey and some preliminary scientific results 
have been discussed only briefly. We expect to use the 
VENGA sample to conduct a series of studies on star for- 
mation, galactic structure assembly, stellar populations, 
chemical evolution, galactic feedback, nuclear activity, 
and the properties of the interstellar medium in massive 
disk galaxies, and to present these results in a series of fu- 
ture publications. As an example of the practical uses of 
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this data, in a recent paper we have presented a study of 
the radial profile in the CO to H2 conversion factor Xco 
and its relation to other parameters like the metallic- 
ity, the gas surface density and the ioniz ation p arameter 
across the disk of NGC 628 (jBlanc et al.|[20Ta) . 

Once all the VENGA data is processed we expect to 
make it publicly available to the community. The rich- 
ness of large IFU datasets like the one we are compiling 
goes beyond the scientific goals of our team. We expect 
VENGA to be a useful resource that will complement the 
wealth of multi- wavelength datasets that the community 
has acquired on nearby spiral galaxies over the last few 
decades. 

The VENGA collaboration acknowledges the gener- 
ous support from the Norman Hackerman Advanced Re- 
search Program (NHARP) ARP-03658-0234-2009, GAB. 
acknowledges the support of Sigma Xi, The Scientific 
Research Society Grant in Aid of Research. NJE 
and ALH were supported by NSF Grant AST-1109116. 
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The VENGA Sample 
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NGC 


3949 


11:53:41.5 


47:51:35 


SA(s)bc 


57 


120 


2.9 x 1.7 


19.1 ± 3.1 


TF 


92 


-22.80 ±0.35 




1 


NGC 


4013 


11:58:31.7 


43:56:48 


SAb 


90 


66 


5.2 x 1.0 


18.9 ± 3.1 


TF 


92 


-23.75 ±0.35 


22.95 


2 


NGC 


4254 


12:18:49.4 


14:25:07 


SA(s)c 


30 


60 c 


5.4 x 4.7 


14.3 ± 3.5 


TF 


70 


-26.07 ±0.15 


21.16 


3 


NCG 


4314 


12:22:32.2 


29:53:47 


SB(rs)a 


28 


145 c 


4.2 x 3.7 


9.7± 3.6 


TF 


87 


-23.81 ±0.15 


21.11 


2 


NGC 


4450 


12:28:29.4 


17:05:05 


SA(s)ab 


43 


175 


5.2 x 3.9 


15.3 ± 2.5 


TF 


74 


-23.87 ±0.35 


21.79 


2 


NGC 


4569 


12:36:50.1 


13:09:48 


SAB(rs)ab 


65 


23 


9.5 x 4.4 


9.9 ±0.2 


STF 


48 


-23.39 ±0.06 


22.10 


1 


NGC 


4826 


12:56:44.3 


21:41:05 


SA(rs)ab 


59 


115 


10.0 x 5.4 


4.4 ±0.1 


TRGB 


21 


-22.87 ±0.03 


20.69 


1 


NGC 


5055 


13:15:49.3 


42:02:06 


SA(rs)bc 


57 


105 


12.6 x 7.2 


9.0 ±0.1 


TRGB 


Id 


-24.16 ±0.03 


21.39 


1 


NGC 


5194 


13:29:53.4 


47:11:48 


SA(s)bc 


20 b 


163 


11.2 x 6.9 


9.1 ± 0.6 


Flow 


44 


-24.30 ±0.15 


21.40 


3 


NGC 


5713 


14:40:11.6 


-00:17:26 


SAB(rs)bc 


28 


10 


2.8 x 2.5 


31.3 ± 2.2 


Flow 


152 


-24.15 ±0.15 


21.36 


1 


NGC 


5981 


15:37:53.4 


59:23:34 


Sc 


90 


110 


2.8 x 0.5 


49.7 ±9.2 


TF 


241 


-24.19 ±0.40 




1 


NGC 


6503 


17:49:27.7 


70:08:41 


SA(s)cd 


71 


123 


7.1 x 2.4 


4.0 ±0.1 


TRGB 


19 


-20.71 ±0.04 


21.08 


2 


NGC 


6946 


20:34:52.3 


60:09:14 


SAB(rs)cd 


32 


60 c 


11.5 x 9.8 


6.1 ± 0.6 


PNLF 


29 


-23.55 ±0.21 


22.93 


2 


NGC 


7479 


23:04:57.1 


12:19:18 


SB(s)c 


42 


25 


4.1 x 3.1 


30.2 ± 5.6 


TF 


146 


-24.20 ±0.40 


22.42 


2 


NGC 


7331 


22:37:04.0 


34:24:56 


SA(s)b 


72 


171 


10.5 x 3.7 


14.5 ± 0.6 


Ceph 


70 


-24.78 ±0.09 


21.51 


1 



a Coordinates, inclination (i), position angle (9), isophotal diameter (cZ 2 5), and effective S-band surface brightness (/jg) taken from RC3 i dc Vaucoulcurs ct al. 1991) except when 
indicated. 

b For NGC5194 we use kinematic inclination angle derived bv lTullvl <1974| >. 

C From Paturel et al. 12000) (NGC1042, NGC3938), Sprinaob et al. 12009) (NGC4254, NGC6946), and .Jarrett et al. (2003) (NGC4314). 

d Distance (L>) methods and references; TRGB: Tip of the red giant branch Uacobs et al.|[2'UU5| ; ITully et al.||2009l except for NGC3351 taken from PTT^zi et alj <200l) ; Ceph: 
Cepheid variables IFreedman et al.| <200l|) , except for NGC2841 taken from rCTaT:ri et al.| <2001| )); TF: HI 21cm Tully-Fisher <Tully et"aTll^QQ3| , for NGC0337 we used the group Tully- 
Fisher distance) except for NGC4252 taken from ^SjiFin gob et all 120091 1; STF: Stellar kinematics Tully-Fisher ICortes eFaTI \2003) ; PNLF: Planetary nebulae luminosity function 
< Herrmann et al Jf^MsT) ; Flow: Derived from redshift, and corrected for peculiar velocities jMould et"ani200o1 , taken from NED). 

6 From PaTrett et all <2003|) except for NGC 1042, NGC 3147, NGC 3949, NGC 5981, NGC 7479, and NGC 7331 taken from PaTrett et alj <2000| ). 
From this work. 
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TABLE 2 
Bulge Structural Parameters 



Object 


B/T 


n Bulg- 


NGC 337 


- 




NGC 628 


0.10 a 


1.35 


NGC 1042 




" 


NGC lUOO 






NGC 2775 


0.61 b 


4.85 


NGC 2841 


0.17 c 


2.97 


NGC zyUo 


0.09 c 


0.42 


NGC o!4/ 


u.zo 


O.OO 


NGC oloo 


0.25 b 


0.56 


NGC olyo 


0.11 a 


5.12 


IN GC OZZ I 






NGC oool 


U.17 


1.51 


IN GC OOZL 


U.lU 


o.zU 


IN v_i V_ OUZ / 


U.UO 


9 QO 
z.yu 


NGC o9oo 


0.07 b 


1.18 


vpn on /in 
NGC o94y 






wpp /I HI Q 
In LtU Wlo 






NGC 4zo4 


n onb 

U.oy 


2.68 


NGC 4ol4 






NGC 4450 


0.17 b 


2.26 


NGC 4569 


0.06 c 


1.90 


NGC 4826 


0.13 c 


3.94 


NGC 5055 


0.26 a 


1.84 


NGC 5194 






NGC 5713 


0.33 b 


1.84 


NGC 5981 






NGC 6503 






NGC 6946 






NGC 7479 


0.09 b 


1.09 


NGC 7331 







a K-band Decomposition, Dong Sz Dc Robcrtis 2006 
b H-band Decomposition, Wcinzirl ct al. 2008 
c V-band Decomposition, Fisher & Drory 2008 
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TABLE 3 

Stellar Masses and Star Formation 
Rates 



Object log(M») log(SFR) Ref. 
M M yr _1 



NGC 


337 


10.2 


0.63 


K03 a 


NGC 


628 


10.3 


0.30 


L09 c 


NGC 


1042 


8.9 


0.15 


T07 b 


NGC 


1068 


10.9 


1.59 


T07 


NGC 


2775 


11.1 


0.06 


T07 


NGC 


2841 


11.1 


-0.70 


K03 


NGC 


2903 


10.7 


0.56 


L09 


NGC 


3147 


11.5 


- 


- 


NGC 


3166 


11.0 


- 


- 


NGC 


3198 


10.4 


-0.07 


K03 


NGC 


3227 


10.8 


- 


- 


NGC 


3351 


10.0 


0.20 


L09 


NGC 


3521 


11.0 


0.38 


L09 


NGC 


3627 


10.7 


0.69 


L09 


NGC 


3938 


10.5 


0.08 


K03 


NGC 


3949 


10.3 






NGC 


4013 


10.7 






NGC 


4254 


11.6 


1.04 


K03 


NGC 


4314 


10.7 


-0.18 


T07 


NGC 


4450 


10.8 


-0.30 


K03 


NGC 


4569 


10.6 


0.28 


K03 


NGC 


4826 


10.4 


0.07 


L09 


NGC 


5055 


10.9 


0.48 


L09 


NGC 


5194 


10.9 


0.88 


L09 


NGC 


5713 


10.9 


0.99 


T07 


NGC 


5981 


10.9 






NGC 


6503 


9.5 


-0.50 


L09 


NGC 


6946 


10.6 


0.96 


L09 


NGC 


7479 


10.9 


1.21 


T07 


NGC 


7331 


11.1 


0.62 


K03 



a K03: Kennicutt ct al. 2003 
b T07: Thilkor ot al. 2007 
c L09: Lee et al. 2009 
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TABLE 4 
VENGA Observing Runs 



Dates Observed Nights Instrumental Setup Observed Galaxies 

08/04/2008 1 red NGC 5194 

11/04/2008 - 11/09/2008 6 red NGC 628, NGC 1068 

01/28/2009 - 01/31/2009 4 red NGC 2903, NGC 3521 

02/01/2009 - 02/03/2009 3 red NGC 1042, NGC 2775, 

NGC 3227, NGC 3949, 
NGC 4314 

03/30/2009 - 04/02/2009 3 red NGC 3351, NGC 4254 

04/17/2009 - 04/19/2009 4 red NGC 4254, NGC 5194 

07/15/2009 - 07/16/2009 2 red NGC 5713, NGC 6503 

07/21/2009 - 07/23/2009 1 red NGC 6503 

09/11/2009 - 09/15/2009 3 red NGC 337, NGC 1068, 

NGC 5981, NGC 6503, 
NGC 7479 

11/09/2009 - 11/15/2009 4 red NGC 628, NGC 1042, 

NGC 1068, NGC 2775 

12/09/2009 - 12/21/2009 12 red NGC 628, NGC 1042, 

NGC 2775, NGC 2841, 
NGC 3166, NGC 3227, 
NGC 3521, NGC 3627 

01/11/2010 - 01/16/2010 4 red NGC 1042, NGC 2841, 

NGC 3147, NGC 3627, 
NGC 4013 

02/14/2010 - 02/18/2010 4 red NGC 1068, NGC 2775, 

NGC 3147, NGC 3198, 
NGC 4013, NGC 4254 

05/18/2010 - 05/20/2010 3 red NGC 3998, NGC 5055 

06/04/2010 - 06/06/2010 2 red NGC 3198, NGC 4450, 

NGC 6964 

07/05/2010 - 07/09/2010 3 red NGC 4450, NGC 6946 

09/01/2010 - 09/07/2010 6 blue NGC 1068, NGC 5981, 

NGC 6503, NGC 6946, 
NGC 7479, NGC 7731 

10/01/2010 - 10/06/2010 5 blue NGC 337, NGC 628, 

NGC 1068, NGC 6946, 
NGC 7731 

11/10/2010 - 11/14/2010 2 blue NGC 628 

12/07/2010 - 12/12/2010 5 blue NGC 628, NGC 1042, 

NGC 2775, NGC 2903 

12/27/2010 - 01/02/2011 5 blue NGC 628, NGC 2775, 

NGC 2841, NGC 3147, 
NGC 3227 

01/27/2011 - 02/02/2011 4 blue NGC 628, NGC 2775, 

NGC 2841, NGC 3147, 
NGC 3166, NGC 3198, 
NGC 3227 

02/07/2011 - 02/10/2011 3 blue NGC 2841, NGC 3147, 

NGC 3166, NGC 3198, 
NGC 3227, NGC 3521 

03/28/2011 - 03/31/2011 4 blue NGC 2775, NGC 3147, 

NGC 3351, NGC 3627, 
NGC 5713 

04/08/2011 - 04/10/2011 1 blue NGC 3351, NGC 3949 

05/04/2011 - 05/09/2011 5 blue NGC 3227, NGC 3949, 

NGC 3938, NGC 5194, 
NGC 5055, NGC 3521, 
NGC 3521, NGC 4826 

06/24/2011 - 07/04/2011 5 blue NGC 4254, NGC 4826, 

NGC 4013 

03/16/2012 - 03/22/2012 4 blue NGC 3521, NGC 4013, 

NGC 4314, NGC 4450 
NGC 4569 
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TABLE 5 

Summary of VENGA Observations of NGC 628 



Pointing Equatorial Coord. Setup Dither Exposure Time N (Seeing) (Transparency) 
a S hours " 



red Dl 4.00 12 2.06 0.87 

PI 01:36:42.45 15:47:04.6 red D2 3.33 10 2.23 0.87 

red D3 3.33 10 2.20 0.89 

blue Dl 3.33 8 2.19 0.71 

blue D2 2.92 7 2.06 0.64 

blue D3 3.33 8 1.55 0.73 



red Dl 3.00 6 1.92 0.65 

P2 01:36:49.45 15:47:04.2 red D2 3.50 7 2.00 0.67 

red D3 3.50 7 1.87 0.68 

blue Dl 5.00 12 1.72 0.68 

blue D2 5.00 12 1.53 0.69 

blue D3 5.00 12 2.23 0.65 



red Dl 8.50 17 2.62 0.63 

P3 01:36:35.51 15:47:05.0 red D2 7.50 15 2.72 0.61 

red D3 7.50 15 2.57 0.68 

blue Dl 5.42 13 2.27 0.49 

blue D2 5.42 13 2.82 0.58 

blue D3 4.58 11 2.18 0.58 
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TABLE 6 
Fitted Emission Lines in NGC 628 



Transition Wavelength Median S/N Fraction - 5<r Fraction - 3<r 

A 



OII] a 


3726.03 


16.4 


0.99 


1.0 


Ollf 


3728.73 


- 


- 


- 


NcIII] 


3868.69 


1.3 


<0.01 


<0.01 


Nelll 


3967.40 


0.7 


<0.01 


<0.01 


H8 




oooy.uu 


9 n 


U. lO 


97 


He 




3970.08 


2.7 


0.18 


0.42 


H<5 




4101.73 


3.0 


0.26 


0.44 


H 7 




4340.47 


8.0 


0.74 


0.89 


[OIII] 


4363.15 


0.9 


<0.01 


<0.01 


Hell 


4685.74 


0.4 


<0.01 


<0.01 


H/9 




4861.32 


23.5 


0.97 


0.99 


[OIII] 


4958.83 


3.4 


0.24 


0.58 


[OIII] 


5006.77 


9.5 


0.89 


0.98 


NI] a 
NI] a 


5197.90 


1.9 


0.01 


0.15 


5200.39 








NII1 


6547.96 


11.3 


0.92 


0.98 


Ha 




6562.80 


47.4 


1.00 


1.00 


Nil 




6583.34 


26.2 


0.99 


1.00 


SII 




6716.31 


16.5 


0.98 


1.00 


SII 




6730.68 


10.4 


0.88 


0.97 



a Since we cannot resolve the [OII]A3727 and [NI]A5200 doublets, we report the 
median S/N and fraction of the observed area in which the lines are significantly 
detected for the sum of the two doublet components. 



